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Abstract. This article is concerned with the derivation and the mathematical 
study of a new mean-field model for the description of interacting electrons 
in crystals with local defects. We work with a reduced Hartree-Fock model, 
obtained from the usual Hartree-Fock model by neglecting the exchange term. 

First, we recall the definition of the self-consistent Fermi sea of the perfect 
crystal, which is obtained as a minimizer of some periodic problem, as was 
shown by Catto, Le Bris and Lions. We also prove some of its properties 
which were not mentioned before. 

Then, we define and study in detail a nonlinear model for the electrons of 
the crystal in the presence of a defect. We use formal analogies between the 
Fermi sea of a perturbed crystal and the Dirac sea in Quantum Electrodynam- 
ics in the presence of an external electrostatic field. The latter was recently 
studied by Hainzl, Lewin, Sere and Solovej, based on ideas from Chaix and 
Iracane. This enables us to define the ground state of the self-consistent Fermi 
sea in the presence of a defect. 

We end the paper by proving that our model is in fact the thermodynamic 
limit of the so-called supercell model, widely used in numerical simulations. 

Describing the electronic state of crystals with local defects is a major issue in 
solid-state physics, materials science and nano-electronics [251 [T71 155] , 

In this article, we develop a theory based on formal analogies between the Fermi 
sea of a perturbed crystal and the polarized Dirac sea in Quantum Electrodynam- 
ics in the presence of an external electrostatic field. Recently, the latter model was 
extensively studied by Hainzl, Lewin, Sere and Solovej in the Hartree-Fock approx- 
imation [TOl HU H31 Q2] , based on ideas from Chaix and Iracane [6] (see also [HQ])- 
This was summarized in the review |14) . Using and adapting these methods, we are 
able to propose a new mathematical approach for the self-consistent description of 
a crystal in the presence of local defects. 

We focus in this article on the reduced Hartree-Fock (rHF) model in which the so- 
called exchange term is neglected. To further simplify the mathematical formulas, 
we do not explicitly take the spin variable into account and we assume that the 
host crystal is cubic with a single atom of charge Z per unit cell. The arguments 
below can be easily extended to the general case. 

In the whole paper, the main object of interest will be the so-called density 
matrix of the electrons. This is a self-adjoint operator < 7 < 1 acting on the 
one-body space L 2 (M 3 ). When 7 has a finite rank, it models a finite number of 
electrons. In the periodic case, the ground state density matrix 7p er has an infinite 
rank (it describes infinitely many electrons) and commutes with the translations 
of the lattice. We will see in the sequel that the ground state density matrix of a 
crystal with a local defect can be written as 7 = 7p Cr + Q, where Q is a compact 
perturbation of the periodic density matrix 7p 0r of the reference perfect crystal. 
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In each of the above three cases (finite number of electrons, perfect crystal, 
defective crystal), the ground state density matrix can be obtained by minimizing 
some nonlinear energy functional depending on a set of admissible density matrices. 
In the case of a crystal with a local defect, the perturbation Q is a minimizer of 
some nonlinear minimization problem set in the whole space R 3 , with a possible 
lack of compactness at infinity. The main unusual feature compared to standard 
variational problems is that Q is a self-adjoint operator of infinite rank. This was 
already the case in \W \ fTf ] fT2 l fT3] . 

The paper is organized as follows. In Section 1, we recall the definition of the 
reduced Hartrcc-Fock model for a finite number of electrons, which serves as a 
basis for the theories of infinitely many electrons in a (possibly perturbed) periodic 
nuclear distribution. Section 2 is devoted to the definition of the model for the 
infinite periodic crystal, following mainly [H [5] (but we provide some additional 
material compared to what was done in [4J [5]). In Section 3, we define a model 
for the crystal with local defects which takes the perfect crystal as reference. In 
Section 4, we prove that this model is the thermodynamic limit of the supercell 
model. 

For the convenience of the reader, we have gathered all the proofs in Section 5. 
Often, the proofs follow the same lines as those in p~0j [Til H21 US] and we shall 
not detail identical arguments. But there are many difficulties associated with the 
particular model under study which do not appear in previous works and which are 
addressed in detail here. 

1. THE REDUCED HARTREE-FOCK MODEL FOR N ELECTRONS 

We start by recalling the definition of the reduced Hartree-Fock model [3JJ for a 
finite number of electrons. Note that the reduced Hartree-Fock model should not 
be confused with the restricted Hartree-Fock model commonly used in numerical 
simulations (see e.g. [S]). We consider a system containing N nonrelativistic quan- 
tum electrons and a set of nuclei having a density of charge p n uc- If for instance 
there are K nuclei of charges z\,...,Zk € N \ {0} located at R\, ...,Rk & R 3 , then 

K 

Pnuc(x) := ^ z k m k (x - R k ), 

k=l 

where mi, ttlk are positive measures on R 3 of total mass one. Point-like nuclei 
would correspond to m k = 5 (the Dirac measure) but for convenience we shall deal 
with smeared nuclei in the sequel, i.e. we assume that for all k = 1...K, m k is 
a smooth nonnegative function such that L 3 m k = 1. The technical difficulties 
arising with point-like nuclei will be dealt with elsewhere. 

The energy of the whole system in the reduced Hartree-Fock model reads [3TJ [5] 

(1.1) ^HF (7) = Tr (-^A^j + X -D { Pl - p mc , Pl - p nuc ) . 

We have chosen a system of units such that H = m = e = j^— = 1 where m and 
e are respectively the mass and the charge of an electron, H is the reduced Planck 
constant and eo is the dielectric permittivity of the vacuum. The first term in 
the right-hand side of (jl.ip is the kinetic energy of the electrons and D(-,-) is the 
classical Coulomb interaction, which reads for / and g in L 6 / 5 (R 3 ) as 

(1.2) D(f, 9 )=[ I iMMd X dy = A,[ 

Jr3J R 3 \x-y\ J R 3 Ife^ 

where / denotes the Fourier transform of /. In this mean-field model, the state of 
the N electrons is described by the one-body density matrix 7, which is an element 
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of the following class 

V N = j 7 € 5(L 2 (M 3 )) | < 7 < 1, Tr( 7 ) = N, Tr (V^A^V^a) < oo 

Here and below, S(Sj) denotes the space of bounded self-adjoint operators acting on 
the Hilbert space $). Also we define Tr((— A) 7 ) := Tr(V— A 7 \/— A) which makes 
sense when 7 6 "P^. The set is the closed convex hull of the set of orthogonal 
projectors of rank N acting on L 2 (R 3 ) and having a finite kinetic energy. Each 
such projector 7 = J^., \ipi){ipi\ is the density matrix of a Hartree-Fock state 

(1.3) * = y>i A • • • AifiN 

in the usual iV-body space of fermionic wavefunctions with finite kinetic energy 

The function p 1 appearing in (11. ip is the density associated with the operator 7, 
defined by p-y(x) — j(x, x) where 7(0:, y) is the kernel of the trace class operator 7. 
Notice that for all 7 S V N , one has p 7 > and yfp^ £ H 1 (M. 3 ), hence the last term 
of (TUT]) is well-defined, since p 1 £ L 1 ^ 3 ) n L 3 (R 3 ) C £ 6 / 5 (R 3 ). 

It can be proved (see the appendix of [3T]) that if N < J^t-Li z k (neutral or 
positively charged systems), the variational problem 

(1-4) J rH F(p n uc,A0=mf{C!(7), TeP"} 

has a minimizer 7 and that the corresponding minimizing density p 7 is unique. 

The Hartree-Fock model [3T] is the variational approximation of the time- independent 
Schrodinger equation obtained by restricting the set of fermionic wavefunctions un- 
der consideration to the subset of functions of the form (|1.3|l . The HF functional 
reads 

the last term being called the exchange energy. As the Hartree-Fock energy func- 
tional is nonconvex, there is little hope to obtain rigorous thermodynamic limits in 
this setting, at least with current state-of-the-art techniques. For this reason, the 
exchange term is often neglected in mathematical studies. 



2. The reduced Hartree-Fock model for a perfect crystal 

In this article, we clamp the nuclei on a periodic lattice, optimizing only over 
the state of the electrons. More precisely we are interested in the change of the 
electronic state of the crystal when a local defect is introduced. To this end, we shall 
rely heavily on the rHF model for the infinite perfect crystal (with no defect) which 
was studied by Catto, Le Bris and Lions in [U[S]. The latter can be obtained as the 
thermodynamical limit of the rHF model for finite systems which was introduced 
in the previous section. This will be explained in Section |4] below. 

Let T = [-1/2, 1/2) 3 be the unit cell. We denote by V* = [-7r,7r) 3 the first 
Brillouin zone of the lattice, and by the translation operator on L 1 2 oc (R 3 ) defined 
by Tku(x) = u(x — k). We then introduce 

Pper = |t e S(L 2 {R 3 )) I < 7 < 1, Vfc e Z 3 , r fc7 = 7 T fe , 

£ Tr L , (r) ((l - A e ) 1 / 2 7e (l - Af)V3) d£ < 00 J 
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where (7§)fer* is the Bloch waves decomposition of 7, see [2"T1 15]: 



7= (27T) 3 

L|(r) = {u e iL( K3 ) I T ku = e-' lk< u, Vfc e z 3 } 

which corresponds to the decomposition in fibers L 2 (R 3 ) = L?(r)d£. For any 
7 G "Ppcr, we denote by ^y^{x,y) the integral kernel of 7^. The density of 7 is then 
the nonnegative Z 3 -periodic function of L 1 1 oc (lR 3 ) n L 3 oc (R 3 ) defined as 

Notice that for any 7 £ Ppcr 

Pl {x)dx = ^3 ^ Tr L 2 (r) ( 7? ) d£, 

i.e. this gives the number of electrons per unit cell. Later we shall add the constraint 
that the system is neutral and restrict to states 7 E "P pC r satisfying 



r 



p~ l {x)dx = Z 



where Z is the total charge of the nuclei in each unit cell. 

We also introduce the Z 3 -periodic Green kernel of the Poisson interaction 
denoted by G\ and uniquely defined by 



j -AGi =4tt( h-1 
G\ Z 3 -periodic, minGi = 

I, K 3 



= 0. 

' R 3 

The Fourier expansion of G\ is 



\k\ 

fc£2n-Z 3 \{0} 1 1 



with c = JpGi > 0. The electrostatic potential associated with a Z 3 -periodic 
density p £ Lj^R 3 ) n L 3 oc (R 3 ) is the Z 3 -periodic function defined as 



(p*r Gi)(x) := J^G 1 {x~y)p{y)dy. 
We also set for any Z 3 -periodic functions / and g 



DaALg) ■= J J Gx{x-y)f{x)g{y)dxdy. 

Throughout this article, we will denote by Xi the characteristic function of the set 
/Cl and by Xi(A) the spectral projector on / of the self-adjoint operator A. 
The periodic density of the nuclei is given by 

(2.1) Mper(aO = ^ Zm(x-R). 



We assume for simplicity that m is a nonnegative function of (R 3 ) with support 
in r, and that J R3 m{x)dx = 1. Hence J r p peT (x)dx = Z, the total charge of the 
nuclei in each unit cell. The periodic rHF energy is then defined for 7 <E Vpcr as 



(2.2) 



£per(7) = ^~)3 j r Tl L 2 ^) [~\ A ^J ^ + \ DGl ^ ~ ^1 Pi ~ ' 
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Introducing 

(2-3) V* cr :-- I - = P,„, ! / ^ ^ Z 

the periodic rHF ground state energy (per unit cell) is given by 
(2.4) 



/ p cr = inf{f p ° er ( 7 ), 7 6 7>p Z CT } 



It was proved by Catto, Le Bris and Lions in [5] that there exists a minimizer 
7per G Ppei *° the minimization problem (|2.4|) , and that all the minimizers of (|2.4p 
share the same density p 7 o . We give in Appendix [5] the proof of the following 

Theorem 1 (Definition of the periodic rHF minimizer). Let Z G N \ {0}. The 
minimization problem (|2.4|) admits a unique minimizer 7 p er . Denoting by 

(2.5) H° eI := - — + (p 7 o M - fi pcr ) * r d, 

f/ie corresponding periodic mean-field Hamiltonian, 7p Cr is solution to the self- 
consistent equation 

( 2 -6) 7pcr = X(-oo,e F ] (H° el ), 

where cf is a Lagrange multiplier called Fermi level, which can be interpreted as a 
chemical potential. 

Additionally, for any tj?gl such that (|2.6p holds, 7 p 0r is the unique minimizer 
on Vper of the energy functional 

7 ^ £per(7) - £F^P 7 - 

Theorem [T] contains three main results that were not present in [5] : first 7p er is 
unique, second it is a projector, and third it satisfies Equation (|2.6p . These three 
properties are crucial for a proper construction of the model for the crystal with a 
defect. 

It can easily be seen that (p 7 o — /i pcr ) *r G\ belongs to L^ oc (IR 3 ). By a result of 
Thomas |34j this implies that the spectrum of H® cr is purely absolutely continuous. 
This is an essential property for the proof of the uniqueness of 7 p 0r . Let (A„(£)) n >i 
denote the nondecreasing sequence of the eigenvalues of (Hp Cr )^. Then 

^p°er) = U A «( r )' ff Per = 7^3 / ^ 

n>l \ n ' •' T * 

The projector 7p er represents the state of the Fermi sea, i.e. of the infinite system 
of all the electrons in the periodic crystal. Of course, it is an infinite rank projector, 
meaning that 

7por = \fk){<Pk\ 
k 

should be interpreted as the one-body matrix of a formal infinite Slater determinant 

* = ipx A ip 2 A • • • A ip k A • • • . 

The fact that 7p 0r is additionally a spectral projector associated with the continuous 
spectrum of an operator leads to the obvious analogy with the Dirac sea which is 
the projector on the negative spectral subspace of the Dirac operator [TO l [TT1 fT2| IT3] . 



Most of our results will hold true for insulators (or semi-conductors) only. When 
necessary, we shall take ZeN\ {0} and make the following assumption: 
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(Al) There is a gap between the Z-th and the (Z + \)-st bands, i.e. < S^+ij 
where and ^z+i are respectively the maximum and the minimum of the Z-th 
ana the [Z + l)-st bands of H° el . 

We emphasize that Assumption (Al) is a condition on the solution 7p Cr of the 
nonlinear problem ()2.4|) . Note that under (Al), one has 7p Cr = X(-oo,e F ] (Hp C r) f° r 
any e F G (£+,£ z+1 ). 

3. The reduced Hartree-Fock model for a crystal with a defect 

In this section, we define the reduced Hartree-Fock model describing the behavior 
of the Fermi sea and possibly of a finite number of bound electrons (or holes) close to 
a local defect. Our model is an obvious transposition of the Bogoliubov-Dirac-Fock 
model which was proposed by Chaix and Iracane [S] to describe the polarized Dirac 
sea (and a finite number of relativistic electrons) in the presence of an external 
potential. Our mathematical definition of the reduced energy functional follows 
mainly ideas from [TU1 [H] . We shall prove in Section 4 that this model can be 
obtained as the thermodynamic limit of the so-called supercell model. An analogous 
result was proved in 13J for the Bogoliubov-Dirac-Fock (BDF) model. 

Assume that the periodic nuclear density /x pcr defined in (|2.ip is replaced by a 
locally perturbed nuclear density ^ pcr + v. The defect v can model a vacancy, an 
interstitial atom, or an impurity, with possible local rearrangement of the neigh- 
boring atoms. The main idea underlying the model is to define a finite energy by 
subtracting the infinite energy of the periodic Fermi sea 7p er defined in the previous 
section, from the infinite energy of the perturbed system under consideration. For 
the BDF model, this was proposed first in [13]. Formally, one obtains for a test 
state 7 

(3.1) £fl + M C!+^Pcr) " = " Tr K°cr(7 - 7°er)) 

"(aOPb-Tgjfo) , , , i f f Ph-i° et ]( x )Ph-i° et ](y) , , 

: : dxdy + - / : : dxdy. 

\x-y\ 2 J M 3 J M 3 \x - y\ 

Of course the two terms in the left-hand side of (|3.1[) are not well-defined because 
/i pcr is periodic and because 7 and 7p er have infinite ranks, but we shall be able 
to give a mathematical meaning to the right-hand side, exploiting the fact that 
Q := 7 — 7p Cr induces a small perturbation of the reference state 7p er . The formal 
computation (|3.ip will be justified by means of thermodynamic limit arguments in 
Section HI 

3.1. Definition of the reduced Hartree-Fock energy of a defect. We now 

define properly the reduced Hartree-Fock energy of the Fermi sea in the presence of 
the defect v. We denote by 6 P the Schatten class of operators Q acting on L 2 (M 3 ) 
having a finite p trace, i.e. such that Tr(|Q| p ) < 00. Note that Si is the space of 
trace-class operators, and that 62 is the space of Hilbert-Schmidt operators. Let 
n be an orthogonal projector on L 2 (R 3 ) such that both n and 1 — n have infinite 
ranks. A self-adjoint compact operator Q is said to be n-trace class (Q S &Y) 
when Q £ & 2 and nQn, (1 - n)Q(l - n) G 61. Its n-trace is then defined as 
Tr n (<2) = Tr(nQn + (1 - n)Q(l - n)). Notice that if Q G 61, then Q G S? for 
any n and Trn(Q) = Tr(Q). See fTUl Section 2.1] for general properties related to 
this definition. In the following, we use the shorthand notation 

Q— ■= 7perQ7per: Q ++ ■= (1 ~ 7^ex)Q(l " 7p er ), 

6? := &1° BI = {Q e 6 2 I Q++ G 61, Q— G 61} and Tr (Q) := Tr 7 o e (Q). 
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We also introduce the Banach space 

Q = {Q £ 6? | Q* = Q, |V|Q e 6 2 , |V|Q++|V| e Si, |V|Q" |V| e Si} , 
endowed with its natural norm 

(3.2) ||Q|| s :=||Q|| e2 + ||0 ++ || 6i + ||0-1| ei 

+ ll|V|0|| S2 + |||V|Q++|V||| 6i + |||V|Q— |V||| Si . 
The convex set on which the energy will be defined is 
(3-3) K := {Q £ Q | - 7 ° cr < Q < 1 - 7p Cr }- 

Notice that /C is the closed convex hull of states Q € Q of the special form Q = 
7 — 7p Cr , 7 being an orthogonal projector on L 2 (R 3 ). Besides, the number Tro(Q) 
can be interpreted as the charge of the system measured with respect to that of the 
unperturbed Fermi sea. It can be proved |10[ Lemma 2] that Tro(Q) is always an 
integer if Q is a Hilbert-Schmidt operator of the special form Q = 7 — 7p Cr , with 7 
an orthogonal projector. Additionally, in this case, Tro(Q) = when ||Q|| < 1. 

Note that the constraint — 7p Cr < Q < 1 — 7p C1 . in (|3.3|) is equivalent [IJ [10] to 
the inequality 

(3.4) Q 2 < Q++ - Q— 

and implies in particular that Q ++ > and Q < for any Q £ K,. 

In order to define properly the energy of Q, we need to associate a density pQ 
with any state Q G JC. We shall see that pq can in fact be defined for any Q e Q. 
This is not obvious a priori since Q does not only contain trace-class operators. 
Additionally we need to check that the last two terms of ()3.ip are well-defined. For 
this purpose, we introduce the so-called Coulomb space 

C := {pe S'QR 3 ) I D{p,p) < 00} 

where D(f,g) — An L 3 \k\~ 2 f(k)g(k)dk was already defined before in (11. 2jl . The 
dual space of C is the Beppo-Levi space C := {V £ L 6 (R 3 ) | £ L 2 (R 3 )} . We 
now use a duality argument to define pq: 

Proposition 1 (Definition of the density pq for Q £ Q). Assume that Q £ Q. 
Then QV £ S? for any V = V 1 + V 2 £ C'+(L 2 (R 3 ) n L°°(E 3 )) and moreover there 
exists a constant C ( independent of Q and V ) such that 

l^o(QV)|<C|Q|| Q (|7i| c , + |V 2 || £a(RS) ). 

Thus the linear form V £ C + (i 2 (]R 3 ) n L°°(M 3 )) ^ Tt (QV) can be continuously 
extended to C + L 2 (IR 3 ) and there exists a uniquely defined function pq £ CnL 2 (R 3 ) 
such that 

W = Vi + V 2 £ C + (L 2 (R 3 ) n L°°(IR 3 )) , ( PQ , Vx) c c , + [ p Q V 2 = Tr {QV). 

The linear map Q £ Q pq £ C n L 2 (R 3 ) is continuous: 

l|Pollc + ll/°ollx-( K 3 ) < C||Q|j Q . 

Eventually when Q £ &i C Sj, then pq(x) — Q(x,x) where Q(x,y) is the integral 
kernel of Q. 

The proof of Proposition [l] is given in Section [5.21 

Assuming that (Al) holds true, we are now in a position to give a rigorous sense 
to the right-hand side of (|3.1[) for 7 — 7°^ = Q £ JC. In the sequel, we use the 
following notation for any Q £ Q: 

(3.5) Tr (iJ p ° cr Q) := Tr(|if p ° cr - K | 1 / 2 (Q++ - Q—)\H° el - nf' 2 ) + kTt (Q) 
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where k is an arbitrary real number in the gap (S^> ^z+i) (this expression will be 
proved to be independent of k, see Corollary [T] below). Then we define the energy 
of any state Q £ K. as 



(3.6) 



S"(Q) := Tr (ff p ° or Q) - D(p Q) u) + \d{ Pq , Pq ). 



The function v is an external density of charge representing the nuclear charge of 
the defect. For the sake of simplicity, we shall assume that v £ L 1 (M 3 )ni 2 (R 3 ) c C 
throughout the paper, although some of our results are true with a weaker assump- 
tion. We shall need the following 

Lemma 1. Assume that (Al) holds true. For any fixed k in the gap (Sj, S^ +1 ), 
there exist two constants C\ , Ci > such that 

(3.7) Cl (l-A)<|#0 er - K |<c 2 (l-A) 

as operators on L 2 (M 3 ). In particular 

- k|V3(i _ A )-V3 <^, \H°- K \-^(1-A)^ <1/Vcl. 



Similarly, (i?p Cr — «;)(1 — A) 1 and its inverse are bounded operators. 

The proof of the above lemma is elementary; it will be given in Section 15.1.11 
By the definition of Q and Lemma [TJ it is clear that the right-hand side of (|3.5[) 
is a well-defined quantity for any Q £ Q and any k £ (S^, T,^ +1 ). By Proposition 
Q] which states that pq £ C for any Q £ Q, we deduce that (|3.6[) is a well-defined 
functional. 

We shall need the following space of more regular operators 
(3.8) Q r := {Q £ Q \ (-A)Q 2 (-A) £ 6 1; (-A)(Q++ - Q )(-A) e 6J 
and the associated convex set 

Kt := K, n Q r . 

The following result will be useful (its proof will be given below in Section [57 



Lemma 2. T7ie space Q T (resp. the convex set JC r ) is dense in Q (resp. in K) for 
the topology of Q. 

Corollary 1. Assume that (Al) holds true. When Q £ Q r , then H® C[ Q £ 
For any Q £ Q, the expression (|3.5p for Tro(-ffp Cr Q) does not depend on k £ 
( s z> s z+i)- IfQ&JC, then 

(3.9) < ciTr((l- A) 1 / 2 Q 2 (1-A) 1 / 2 ) 

< Cl Tr((l - A) 1 / 2 (Q++ - Q— )(1 - A) 1 / 2 ) 

< Tr (i/ p or Q)- K Tr (Q) 

< c 2 Tr((l - A) 1 / 2 (Q++ - Q— )(1 - A) 1 / 2 ) 
where C\ and Ci are given by Lemma\l\ 

Proof. Let Q £ Q t and k £ (S+ Then ((if° er -«;)Q)++ = \H$ el - K \Q ++ = 

\H° el - - A)-!(l - A)Q++ £ 6i by Lemma □ and the definition of Q r . A 
similar argument for ((-ffp er — k)Q)~~ proves that H® eT Q £ 6°. Then for any 
Q £ Qr, (13. 9|) is a straightforward consequence of ()3.7() and (|3.4|) . We conclude 
using the density of Q r in Q and the density of ZC r in /C. □ 



The following is an adaptation of [TU1 Thm 1]: 
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Corollary 2. Let v e L 1 {R 3 ) n L 2 (R 3 ), 2eN\{0} and assume that (Al) holds. 
For any k G (S^> ^z+i)> orie ^ as / or some di, <fe > 

(3.10) VQe/C, f(Q) - kTto(Q) > <*i ( |Q ++ || 6l + || Sl 

-+IV7III _l_lllV7ln IV7lll l.L/7-Anll 2 _LlllY7lnll 2 ^ ' 



+ |||V|Q ++ |V||| ei + |||V|Q— |V||| Si J +d 2 [\\Qf 62 + |||V|Qir 62 j - -D(y,u) 

Hence E v — ftTro is bounded from below and coercive on IC. Additionally, when 
v = 0, Q Ht £ (Q) — K>Ttq(Q) is nonnegative, being its unigue minimizer. 

Proof. Inequality (|3.10p is a straightforward consequence of (13. 9|) and the fact that 
£>(•,•) defines a scalar product on C. The rest of the proof is obvious. □ 

Remark 1. The energy £ U (Q) measures the energy of a state 7 = 7p 0r + Q with 
respect to that of 7p er . Thus the last statement of Corollary [2] is another way of 
expressing the fact that 7p er is the state of lowest energy of the periodic system 
when there is no defect. 

3.2. Existence of minimizers with a chemical potential. In view of Corol- 
lary [21 it is natural to introduce the following minimization problem 



(3.11) E" := inf{£"(Q) - e F Tr (Q), Q G K} > -00 



for any Fermi level €f G C^z^z+i)- ^ nc following result is proved in Section 157 
following ideas from |llj : 

Theorem 2 (Existence of minimizers with a chemical potential). Let v E L 1 (M 3 )P\ 
L 2 (R 3 ), Z G N\{0} and assume that (Al) holds. Then for any ep £ (E^, S^+i)' 
there exists a minimizer Q G K, for l|3.11[) . Problem (13 . may have several min- 
imizers, but they all share the same density p — pq. Any minimizer Q of (|3.1ip 
satisfies the self- consistent eguation 

[ ., v y ) j Q = X(-oo,e F )(H Q ) - ~/° cr + 5, 



1 C„ - _ ,A * I . 1-1 



* 



where 8 is a finite rank self-adjoint operator satisfying < S < 1 and Ran(<5) C 
ker (^Q - 

Remark 2. It is easily seen that (pq — v) * \ ■ | _1 is a compact perturbation of 
H° eT , implying that Hq is self-adjoint on P(i?° or ) = T>(-A) = H 2 {R 3 ) and that 
a ess (HQ) = cr(Hp el ). Thus the discrete spectrum of Hq is composed of isolated 
eigenvalues of finite multiplicity possibly accumulating at the ends of the bands. 

Recall that the charge of the minimizing state Q obtained in Theorem[2]is defined 
as Tr (Q). Similarly to [TD], [TTJ , it can be proved by perturbation theory that for 
any fixed ep, there exists a constant C(ep) such that when D(y,v) < C(ep), one 
has ker(77g — ep) — {0} and Tiq(Q) = 0, i.e. the minimizer of the energy with 
chemical potential ep is a neutral perturbation of the periodic Fermi sea. 

For a fixed external density v and an adequately chosen chemical potential ep, 
one can have Tro(Q) ^ meaning either that electron-hole pairs have been cre- 
ated from the Fermi sea, and/or that the system of lowest energy contains a finite 
number of bound electrons or holes close to the defect. In the applications, one 
will usually have for a positively charged nuclear defect (y > 0) that the spectrum 
of Hq contains a sequence of eigenvalues converging to the bottom £^+1 °f t ne 
lowest unfilled band (conduction band), and that e F is chosen such that exactly q 
eigenvalues are filled, corresponding to q bound electrons: 

(3.13) Q= (x(-ocX){H Q ) ~ j° cr ) + (xp,e F )(HQ) + 5) :=Q P oi+7c-, 
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pol 



Q 

Figure 1. Decomposition Q — Q po \ + -f e - for not too strong a 
positively charged nuclear defect (v> 0). 

where we have chosen as a reference the center of the gap 



2 

For not too strong a defect density v, one has ker(77g — E) = {0} and Tr (Q po i) = 0. 
Hence Tr(7 e -) = g. Let us assume for simplicity that 6 = and that geN \ {0}. 
Then 
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(3-14) 7c- = XV,e F ){H Q ) = ^ \<Pn){<p n \ 

n=l 

where (<p n ) are eigenfunctions of ffg corresponding to its g first eigenvalues in 

(3.15) #Q<Pn = A„(y9 n . 

Notice that 

(3.16) H Q = -A/2 + (p 7c _ - !/) * | • |- x + V po , 
where 

^pol = (P 7 or - Mpor) *r Gi + p Qpo] * I • I" 1 
is the polarization potential created by the self-consistent Fermi sea and seen by the 
q electrons. Thus the q electrons solve a usual reduced Hartree-Fock equation (|3.15p 
in which the mean-field operator (|3.16p additionally contains the self-consistent 
polarization of the medium. 

The interpretation given in the previous paragraph is different if the positive 
density of charge v of the defect is strong enough to create an electron-hole pair 
from the Fermi sea. 

We end this section by specifying the regularity of solutions of (|3.12|) . The proof 
is given in Section 15.41 

Proposition 2. Let v € L^M 3 ) n L 2 (R 3 ), Z e N \ {0} and assume that (Al) 
holds. Any Q € K, solution of the self- consistent equation (|3.12[) belongs to IC r . 

Remark 3. Notice that it is natural to wonder whether Q e &i, which would in 
particular imply that pq G L 1 (]R 3 ). This is known to be false for the Bogoliubov- 
Dirac-Fock model studied in [101 HH [12 [13] . We do not answer this question for 
our model in the present paper. 

3.3. Existence of minimizers under a charge constraint. In the previous 
section, we stated the existence of minimizers for any chemical potential in the 
gap of the periodic operator -ffp Cr , but of course the total charge Tro(Q) of the 
obtained solution was unknown a priori. Here we tackle the more subtle problem 
of minimizing the energy while imposing a charge constraint. Mathematically this 
is more difficult because although the energy £ U (Q) is convex on K, and weakly 
lower semi-continuous (wise) for the weak-* topology of Q (as will be shown in the 
proof of Theorem [2]), the 7p er -trace functional Q E JC i— > Tro(Q) is continuous but 
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not wise for the weak-* topology of Q: in principle it is possible that a (positive 
or negative) part of the charge of a minimizing sequence for the charge-constrained 
minimization problem escapes to infinity, leaving at the limit a state of a different 
(lower or higher) charge. In fact, we can prove that a minimizer exists under a 
charge constraint, if and only if some binding conditions hold, the role of which 
being to prevent the lack of compactness. 

As explained above, imposing Tro(Q) = q should intuitively lead (for a suffi- 
ciently weak defect density v) to a system of q electrons coupled to a polarized 
Fermi sea. Notice that we do not impose that q = J K3 z/, i.e. our model allows a 
priori to treat defects with non-zero total charge. 

As usual in reduced Hartree-Fock theories, we consider the case of a real charge 
constraint q £ K: 



(3.17) E"{q) := inf{£"(Q), Q £ AC, Tr n (Q) = q}. 



When no defect is present, E (q) can be computed explicitly: 

Proposition 3 (Defect-free charge-constrained energy). Let Z e N \ {0} and as- 
sume that (Al) holds. Then one has 

F o (n \ = / S z+i? when g> 
W \ Z+q when q < 0. 

The minimization problem (|3.17|) with v = has no solution except when q = 0. 

The proof of Proposition [3] is given in Section f5. 61 We now state the main result 
of this section, which is directly inspired from |12j : 

Theorem 3 (Existence of minimizers under a charge constraint). Let v £ L 1 (M 3 )n 
L 2 (M 3 ), Z £ N\ {0} and assume that (Al) holds. The following assertions are 
equivalent: 

(a) Problem (|3.17p admits a minimizer Q; 

(b) Every minimizing sequence for (|3.17|) is precompact in Q and converges 
towards a minimizer Q of (|3.17p . - 

(c) V<z'eK\{0}, E v (q) <E tJ (q-q')+E°(q'). 

Assume that the equivalent conditions (a), (b) and (c) above are fulfilled. Ln 
this case, the minimizer Q is not necessarily unique, but all the minimizers share 
the same density p — pq. Besides, there exists tp € [E^j^z+iI su °h that the 
obtained minimizer Q is a global minimizer for E" p defined in p. lip . Lt solves 
Equation (|3.12p for some < 8 < 1 with Ran(<5) C ker(iJg — ep). The operator 5 
is finite rank if cf € (^z'^z+i) an< ^ trace-class if ep £ {T,'^,T,^ Z+1 }. 

Additionally the set of q's in K satisfying the above equivalent conditions is a 
non-empty closed interval ICR. This is the largest interval on which q <— > E v (q) 
is strictly convex. 

Remark 4. One has I = {Tr (<3), Q min. of E^, e F £ [S+,E Z+1 ]}. Hence I ^ 
by Theorem [U 

Theorem [3] is proved in Section 15.81 Many of the above statements are very 
common in reduced Hartree-Fock theories and not all the details will be given (see, 
e.g. jH]). The difficult part is the proof that (b) is equivalent to (c), for which we 
use ideas from [12"] , 

Conditions like (c) appear classically when analyzing the compactness proper- 
ties of minimizing sequences, for instance by using the concentration-compactness 
principle of P.-L. Lions [23]. They are also very classical for linear models in which 
the bottom of the essential spectrum has the form of the minimum with respect 
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to q' of the right hand side of (c) , as expressed by the HVZ Theorem [THl [351 . 
Assume for simplicity that q > and that Q can be written as in (|3.13|) and (13.14|) . 
When < q 1 < q, (c) means that it is not favorable to let q' electrons escape to 
infinity while keeping q — q' electrons near the defect. When q' < 0, it means that 
it is not favorable to let \q'\ holes escape to infinity, while keeping q+ \q'\ electrons 
near the defect. When q' > q, it means that it is not favorable to let q' electrons 
escape to infinity, while keeping q' — q holes near the defect. 

In this article we do not show when (c) holds true. Proving (c) usually requires 
some decay property of the density of charge pq for a solution Q of the nonlinear 
equation (|3.12|) . In particular, knowing that pq £ L 1 (M 3 ) would be very useful (see 
Remark 13]). We plan to investigate more closely the decay properties of pq and the 
validity of (c) in the near future. At present, the validity of a condition similar 
to (c) for the Bogoliubov-Dirac-Fock model was only proved in the nonrelativistic 
limit or in the weak coupling limit, see |12j . 

4. Thermodynamic limit of the supercell model 

As mentioned before, we shall now justify the model of the previous section by 
proving that it is the thermodynamic limit of the supercell model. 

Let us emphasize that there are several ways of performing thermodynamic lim- 
its. In [5], the authors consider a box of size L, Al := [— L/2, L/2) 3 , and assume 
that the nuclei are located on Z 3 n A^. Then they consider the rHF model of Section 
Q]for TV electrons living in the whole space, with N = ZL 3 chosen to impose neu- 
trality. Denoting by p^ the ground state electronic density of the latter problem, 
it is proved in Thm 2.2] that the energy per unit cell converges to I^ eT , and that 
the following holds: 

y feez 3 nA L 

weakly in Hl oc (M. 3 ), strongly in i{ ) oc (R 3 ) for all 2 < p < 6 and almost everywhere 
on K 3 when L — > oo. Let us recall that 7p er and 7p er are defined in Section^ 

Another way for performing thermodynamic limits is to confine the nuclei and 
the electrons in a domain with ^ 00 > by means of Dirichlct boundary 
conditions for the electrons. The latter approach was chosen for the Schrodinger 
model with quantum nuclei in the canonical and grand canonical ensembles [35] by 
Lieb and Lebowitz in the seminal paper [TU] (see also [H]), where the existence of a 
limit for the energy per unit volume is proved. The crystal case in the Schrodinger 
model was tackled by Fefferman [5] in the same spirit. We do not know whether 
Fefferman's proof can be adapted to treat the Hartree-Fock case. 

Another possibility, perhaps less satisfactory from a physical viewpoint but more 
directly related to practical calculations (see e.g. [8J), is to take Hl = Az and 
to impose periodic boundary conditions on the box A^. Usually the Coulomb 
interaction is also replaced by a (LZ 3 )-periodic Coulomb potential, leading to the 
so-called supercell model which will be described in detail below. This approach has 
the advantage of respecting the symmetry of the system in the crystal case. It was 
used by Hainzl, Lewin and Solovej in [T2] to justify the Hartree-Fock approximation 
of no-photon Quantum Electrodynamics. The supercell limit of a linear model for 
photonic crystals is studied in [3"2"] . 

Of course the conjecture is that the final results (the energy per unit cell and 
the ground state density of the crystal) should not depend on the chosen thermo- 
dynamic limit procedure. This is actually the case for the reduced Hartree-Fock 
model of the crystal. See [15] for a result in this direction for a model with quantum 
nuclei. 
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Let us now describe the supercell model. For L £ N \ {0}, we introduce the 
supercell = [— L/2,L/2) 3 and the Hilbert space 

L 2 per (A L ) = S Lf oc (R 3 ) | <p (LZ 3 )-periodic} . 

We also introduce the LZ 3 -periodic Coulomb potential Gl defined as the unique 
solution to 



AG, 



47T 



Gl L Z 3 -periodic, min R 3 Gl =0 
It is easy to check that Gl(x) — L~ 1 Gi(x/L) and that 

c v-^ 4f 1 



ik-x 



For any LZ 3 -periodic function g, we define 



(g*A L G L )(x) := G L (x -y)g(y)dy, 

Dc L {f,g)-= / G L {x-y)f{x)g{y)dxdy. 
Ja l Ja l 

An admissible electronic state is then described by a one-body density matrix 7 in 

'Pscx = {7 G 6i(4 er (Ai)) I 7* =7, < 7 < 1, Tr L 2 cr(Ai) (-A 7 ) < +00J. 

Any 7 G Psc.l has a well-defined density of charge p-y(x) — j(x,x) where 7(2,?/) 
is the kernel of the operator 7. Notice that 7(2 + Lz, y + Lz') = j(x,y) for any 
z,z' £ Z 3 , which implies that p 1 is LZ 3 -periodic. 

Throughout this section, we use the subscript 'sc' to indicate that we consider 
the thermodynamic limit of the supercell model. 

4.1. Thermodynamic limit without defect. Because our model with defect 
uses the defect-free density matrix of the Fermi sea as a reference, we need to start 
with the study of the thermodynamic limit without defect. We are going to prove 
for the supercell model a result analogous to [5J Thm 2.2]. 

The reduced Hartree-Fock energy functional of the supercell model is defined for 
7 G V sc ,l as 



£ s °cL(7) = Tr iMAi) (-IaJ)+Idc 



where we recall that p, pel (x) = J2rez 3 Zm{x — R) is a Z 3 - (thus LZ 3 -) periodic 
function. The reduced Hartree-Fock ground state energy for a neutral system in 
the box of size L is then given by 



(4.2) 



J S ° C . L = illf (^(7), 7 G V SC ,L, [ P 1 = [ Mper = ZU 
I JA L JA L 



Let us recall that L® er , 7p Cr 



and i?p 0r are defined in Section^ In Section we 



prove the following 

Theorem 4 (Thermodynamic limit of the defect-free supercell model). Let Z > 0. 
i) For all L £ N\{0}, the minimizing problem I® c L has at least one minimizer, and 
all the minimizer s share the same density. This density is Z 3 -periodic. Besides, 
there is one minimizer 7° c L of (|4.2[) which commutes with the translations t^, 
k £ Z 3 . 
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ii) The following thermodynamic limit properties hold true: 

• (Convergence of the energy per unit cell). 

l Q 

lim „ J pcr' 

• (Convergence of the density). 

(4.3) v^cT ~" weakl y in ff ioc( R3 )> 

p 7 o ^ — > p 7 o strongly in L|' oc (R 3 ) for 1 < p < 3 and a.e.; 

• (Convergence of the mean-field Hamiltonian and its spectrum). Let 

seen as an operator acting on L 2 (R 3 ). Then, for all L S N\ {0}, _ff° c L — i?p Cr is a 
bounded operator and 

i:™ II zltO 

,L per | 



L — >oo 



lim II M r -if° || =0. 



Denoting by (A„(£))„ e N\{0} ^ e nondecreasing sequence of eigenvalues of (H® c 
for ^ eT*, one has 

(4.4) lim sup sup |A£(0-An(0| =0 

- L ^ oo riGN\{0}5Gr* 

where (A n (£))„>i are i/ie eigenvalues of (H® eI )£ introduced in Theorem]]] 

Hi) Assume in addition that (Al) holds. Fix some €f € (S^S^+i)- Then for L 
large enough, the minimizer 7g C L of I® c L is unique. It is also the unique minimizer 
of the following problem 

(4-5) I° sc ,L,e P ■■= mf {4°c,l(7) " ^i^^ L ){i), 7 G 7> SC , L } . 

Remark 5. Notice that some of the above assertions are more precise for the 
supercell model than for the thermodynamic limit procedure considered in Thm 
2.2] (compare for instance (|4.3[) with l|4.ip ). This is because the supercell model 
respects the symmetry of the system, allowing in particular to have a minimizer 
7° c L in the box of size L 3 which is periodic for the lattice Z 3 . For an insulator, the 
uniqueness of j® c L for large L and the convergence properties of Hi ) are also very 
interesting for computational purposes. 

4.2. Thermodynamic limit with defect. We end this section by considering 
the thermodynamic limit of the supercell model with a defect. Recall that v € 
L^R 3 ) n L 2 (R 3 ) C C is the density of charge of the defect. First we need to 
periodize this function with respect to the large box Al, for instance by defining 

vl{x) := 5^(l Ai i/)(--L*). 

The reduced Hartree-Fock energy functional of the supercell model with defect is 
then defined for 7 £ V sc ,l as 



£Tc,l(7) = Tr L2 cr (A L ) ( J + 2 Dgl ^ ~ ^P cr ~ Vh ' p i ~ M P° r ~ Vjj } ' 



For ef e {^zi ^z+i)> we consider the following minimization problem 
(4.6) 



!s C ,L,e P = inI {^sc,l(i) - £ F Tr I?„(Ai)(7), 7 G P S c,£ j 



We recall that 7p er is defined in Section that E% and Q are defined in Section 
and that I° c L ep is defined in Section 14.11 In Section 15.101 we prove the 
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Theorem 5 (Thermodynamic limit of the supercell model with defect). Let Z € 
N\ {0}. Assume that (Al) holds and fix some ep € Then one has 

(4.7) lim (j^ £iCF - I° C ^ F ) = E: f - I v (W - Mpor ) * r G x ) + \d{v, v). 

Additionally, if 7^ c L denotes a minimizer for (|4.6p . rten one has, up to extrac- 
tion of a subsequence, 

(lsc,L ~ 1sc,l){x, y) -> QO, y) 

weakly in Hl oc (JS? x R 3 ) and strongly in L^ oc (R 3 xl 3 ), where Q is a minimizer of 
(|3.11[) , as obtained in Theorem^ Besides, 

weakly in Lf oc (M. 3 ), where p is the common density of all the minimizers of (|3.1ip . 

Remark 6. In numerical simulations, the right-hand side of (|4.7[) is approximated 
^ lee he — ^sc l e for a given value of L. This approach has several drawbacks. 
First, the values of L that lead to tractable numerical simulations arc in many 
cases much too small to obtain a correct estimation of the limit L — > oo. Second, 
it is not easy to extend this method for computing E" , to the direct evaluation of 
E u (q) for a given q (i.e. the energy of a defect with a prescribed total charge). The 
formalism introduced in the present article (problems (|3.1ip and (|3.17p ) suggests 
an alternative way for computing energies of defects in crystalline materials. A 
work in this direction was already started [3] . 



5. Proof of the main results 

Unless otherwise stated, the operators used in the following proofs are considered 
as operators on L 2< 



5.1. Useful estimates. We gather in this section some results which we shall need 
throughout the proofs. We start with the 

5.1.1. Proof of Lemma[l[ Recall that -ffp Cr = —A/2 + V pf}r with V pcr := (p 7 o — 
AVr) *r Gi € L°°(l 3 ). Thus \H° el - k\ > H° el - K > -A/2 - C for some large 
enough constant C . On the other hand, as k € (£j, , 1 ), there exists a > such 
that |-ffp er — «| > a. This implies that 

\H° er -k\> max(- A/2 - C, a) > c x (1 - A) 

for some constant c\ > 0. The proof of the upper bound in ()3.7p is straightforward. 

Then (-A/2 + c) -1 - k + c) = 1 + (-A/2 + c) - ^^ - «) is a bounded 
invertible operator for c large enough, since 

IK-A/2+ C )- i (vw-'0ll < iVri£ ; + |K| . 

Thus (i?p Cr — k + c) _1 (— A/2 + c) is bounded for a well-chosen c ^> 1, which clearly 
implies that 

« r - n)-\-A + 1) = g Per -^ + C _ k + c) -l ( _ A/2 + c) . "A + i 



ff° er -K ^ por ■ > v ' ' '-A/2 + c 
is also bounded, together with its inverse. □ 
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5.1.2. Some commutator estimates. Throughout this paper, we shall use Cauchy's 
formula to express the projector 7p er : 



(5.1) 



T P ° er =^- / (z - H^)- 1 dz, 



where ^ is a fixed regular bounded closed contour enclosing the lowest Z bands of 
the spectrum of -ffp Cr . 

The following result will be a useful tool to replace the resolvent (z — _ffp Cr ) _1 
with the operator (—A + 1) _1 which will be easier to manipulate. Its proof is the 
same as the one of Lemma ffl 

Lemma 3. The operator B(z) := (z — Hp eI ) (— A + l) and its inverse are bounded 
uniformly in z G ^€ . 

The next result provides some useful properties of commutators: 

Lemma 4. The operators [jp eT , A] and (1 — A) [7p er , | V|] (1 — A) are bounded. 

Proof. The boundedness of [7p er , A] follows from (|5.ip and from the fact that [(z — 
_ffp Cr ) _1 , A] is bounded uniformly in z £ ^ by Lemma [31 

Using again fSTT]) . it suffices to prove that (1 - A)[(z - H^y 1 , |V|](1 - A) is 
bounded uniformly in z G ^ to infer that (1 - A) [ 7 ° cr , |V|] (1 - A) is bounded. In 
order to prove the uniform boundedness of (1 — A)[(z — iJp er ) _1 , |V|](1 — A), we 
use the formal equality 

(5.2) [(z - A)-\B] = (z- A)- X [A, B](z A)~\ 
We thus obtain 

(5.3) (1 - A)[(z - H° pa r\ |V|](1 - A) = B(z)* [M,V per ] B(z). 
Using (|5.ip and lemma [3 we obtain 

||(l-A)[ 7 ° cr) |V|](l-A)|| <C\\[\V\,V pei ]\\<C\\VV vei \\ Laom . 

□ 

Lemma 5. Let V = Vi + V 2 with Vi G C and V 2 G i 2 (M 3 ). Then [y° el , V] G ©2 
and there exists a positive real constant C such that 

||[7per,V]|| 6a <C(||7l|| C , + ||V a L a(R 3 ) ). 

Proof. Formulas (|5Tj) and (fO]) lead to 

= -^f^B{z){-A + l)- 1 [^V 2 ]{-A + l)- 1 B{zrdz 
= ~7T / B(z)((-A + l)- 1 A)V 2 (-A + l)- 1 B(zydz 

+ -?- / S(z)(-A + l)- 1 U 2 (A(-A + l)- 1 ) J B(z)*dz. 

As (— A+l) -1 and (— A + 1) _1 A are bounded operators, we obtain, using Lemma[21 
||[7p c r ,V2]|| e2 <C||(-A + l)-V 2 || S2 , 
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for some constant C independent of V% . Likewise, 

= - 1 h L B{z) {{ - A + iyldxi) £ ( - A + i)_lflw * dz 



which implies 



I [-ySer. vi] II - < c II (- A + l)" 1 VVi I 



'per' -LJ II e 2 - H v 7 x II © 2 

for some constant C independent of V\. We then use the Kato-Seiler-Simon in- 
equality (see [20J and [301 Thm 4.1]) 

(5.4) Vp > 2, \\f(-iV)g(x)\\ 6p < {2*)-^ \\g\\ LP(R3) |/| £P(R8) 

to infer 

(5-5) ||[^er,^]|| 6a <C"|V 2 || ia(R , )) 

(5.6) Hbper^ilIU < c"||VFi|| i2(a 3 ) = c"|y 1 i c ,. 

□ 

Lemma 6. The operator [ \H® el — k\, |V| ] is bounded for any k in the gap (S^, ^z+i)- 

Proof. We have \H° pm . - k\ = - «) 7 ° er + - «)(1 - 7° cr ) and thus 

[|^er-«l,|V|] = -2(ff° er -K)[7° er ,|V|] + [|V|,y per ](27° er -l) 

= 2( J B( K )*)- 1 (1 - A)[ 7p ° cr , |V|] + [|V|, y por ](2 7 ° cr - 1) 

which gives the result since ||[|V|, Vp Cr ]|| < || Wperlio^RS) and (1 - A)[7° cr , |V|] is 
bounded by Lemma [U □ 

5.2. Proof of Proposition [U Let V = Vi + V 2 where Vi G C and V 2 G L 2 (M 3 ) n 
L°°(R 3 ), and Q G Q. Notice that 

(QV) ++ = Q ++ V(y%) x + Q + -{ 7 peI ,V}( 1 ° peI ) x , 

(5-7) (QV) = Q—V^ ex - Q- + [7p Cr ,^7pcr, 

where (7p er ) = 1 — 7p Cr . We only treat the (QV) term, the argument being the 
same for (QV) ++ . 

First we write Q — Vj° eI = Q — (1 + |V|)(1 + |V|) _1 F7° er and notice that 
(1 + |V|) _1 y is bounded since V% G L°°(R 3 ) by assumption and 

||(1 + IVI)- 1 ^ <CIV4 LB < C llWrL, = C \Vx\c, 

by the Kato-Seiler-Simon inequality (|5.4j) and the critical Sobolev embedding of 
iJ 1 (R 3 ) in L 6 (R 3 ). This proves that Q 1^7p 0r is a trace-class operator. Thus the 
following is true: 

|Tr(Q-y 7p ° or )| = |Tr(Q V)| 

= |Tr((l + |V])Q— (1 + |V|)(1 + IVD-Va + IVI)- 1 )! 

< HQl^Ki + lvD-Mi + lvi)- 1 !^. 
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Then 

IKi+ivir^a+ivir 1 !^ < ki+ivd-ViH^iki+ivd- 

< C[Vi\ L e <C|Vi|| c ,, 

and 

Ki + ivir^a+ivir 1 !^ < la + ivir^Ki + ivi)- 1 !^ 

2 



< 



(i + |vmv 2 | 1/2 _ <ci^| i2 

e 4 



Hence, 



|Tr(Q— V 7 ° )| < C IQIU (|ViU + ||V 2 || La )- 



y pe jl ^ hwuq y vi lie ' l|f2ii L 2j 
For the second term of (15. 7p . we just use Lemma[5]which tells us that Q h [7p Cr : ^l7per £ 
6i since Q~ + G 6 2 and [7p Cr , V] G 6 2 . Additionally 

|Tr(0- + [7p° CI .,^ 7p CI .)| < ||Q- + [7p CI .,^7p cr || ei < C||Q- + | S2 (||Vi|| c , + ||y 2 || i2 ). 
The end of the proof of Proposition [T] is then obvious. □ 

5.3. Proof of Lemma [2j Let Q G Q. For e > 0, we introduce the following 
regularization operator 

(5.8) + 
and set 

Qe '■= ReQRe- 

Notice first that Q e G Q r - Indeed, using the same notation as in Lemma [3 we 
obtain 

fl 1 Hr^af E Hr^af E 

(1 - A)i? £ = (1 - A)(H° - S)- 1 £2= _ ^ r 

which shows that ||(l-A)i? £ || < e- 1 ||B(E)*||. Similarly, ||i? £ (l-A)|| < £- 1 ||B(E)||. 
As R e commutes with 7 p er , we infer 

(1 - A)Q e — (1 - A) = (1 - A)R e Q—R e (l - A) G 6 X . 

Likewise, (1 - A)Q++(1 - A) G &\. Then we show that Q e G K T C K. when 
Q e K.. To prove this, we use the fact that — 7° < Q < 1 — 7 p er is equivalent to 
Q 2 < Q++ - Q— (see Section GDJ . As ||i? £ || < 1, we obtain 

(5.9) (Q £ ) 2 = R t Q{R t ) 2 QR t < R t Q 2 R e < R e (Q++-Q—)R e = (Qe) ++ -(Qe^ 

where we have used that {Re) 2 < 1 and that 7p 0r commutes with R e . Hence, it only 
remains to prove that Q t — > Q for the Q-topology as e — > 0, for any fixed Q G Q- 
We shall need the 



Lemma 7. For any 1 < p < 00 and any fixed Q G & p , one has 
(5.10) 

Proof. Notice that 



(5.10) lim\\R e Q - Q\\ 6p =0 



e|i?<L - El 
iL - 1 = — pG1 



1 + - E| 

satisfies \\R e — 1|| < 1. Hence \\(R e — l)Ql 6p < ||Q|| e . By linearity and density of 
"smooth" finite rank operators in 6 P for any 1 < p < 00, it suffices to prove (|5.10[) 
for Q = I/) (/| with / G i? 2 (K 3 ). Then 

BC-Re - l)l/X/ll ei < II (-Re ~ 1)/|| L2 l/lb < e |K cr - E|/|| L2 ||/|| L2 

which converges to as e — > and controls ||(i? £ — l)|/)(/|||g for 1 < p < 00. □ 
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We are now able to complete the proof of Lemma [2J First, by ([57 
[R e ,\V\] = -eR e [|ff p ° cr -£|,|V|] R e 
and therefore by Lemma [6] there exists a constant C > such that 
(5.11) || [i? e ,|V|] || <Ce. 

Hence, lim c ^o II [Re, |V|] || = 0. Compute now for instance 

\V\(ReQRe-Q)~~M = \V\((R e -l)Q—Re + Q — (Re-l))M 

= [\V\,Re]Q—[Re, |V|] + [\V\,R e }Q— \V\R e 

+(R e - 1)|V|Q— [R e , |V|] + \V\Q—[R e , |V|] 
+(R, - 1)|V|Q — |V|i? e + |V|Q— |V|(-R« - 1). 

Applying either (|5.11[) or Lemma [7] to each term of the previous expression allows 
to conclude that 

lim|||V|(Q--Q--)|V||| =0. 
The proof is the same for the other terms in the definition of ||-||g. Q 

5.4. Proof of Proposition [2j regularity of solutions. Let Q £ Q be of the 

form 

Q = X(-oo, eF )(#p Cr +V)- 7pcr + « 

where < 6 < 1 is a finite rank operator with Ran(i5) C ker(ffp 0r + V — ep) and 
V = p * | • I" 1 for some p £ L X {M?) n L 2 (M 3 ) (in our case p = p Q - v). Note that 
V A eC'nL°°(M 3 ) since 

I^L- < (2-r 3/2 / \v(k)\dk = c f Mi* 



n 2 ; vy R 3 ifcpa + ifcp), 

Since ker(if£ er + 1/ - e F ) C P(if° er + V) = = ff 2 (M 3 ), it is clear that 

the finite rank operator 8 satisfies (1 — A)<5(1 - A) € Si. Thus, up to a change of 
ep, we can assume that ker(Hp C1 . + V — ep) = {0} and that 6 = 0: 

Q-X(-oo,e F )(^per + ^)-7per- 

Then we remark that Q 2 = Q ++ —Q , hence we only have to prove that (1— A)Q £ 

e 2 . 

Let ^ be a smooth curve enclosing the whole spectrum of H® el . + V below ep. 
Since V £ L°°(R 3 ) and \Hp Cr + V — z\ > c > uniformly in z £ we can mimic 
the proof of Lemma [3] and find that 

(5.12) sup ||(1- A)(JZ° . + V-z)- 1 \\ < oo. 



We then use Cauchy's formula (|5.ip and iterate the resolvent formula 

(z - H° pcl . - V)- 1 = {z- H^)- 1 + (z- H° pcl . - V)- l V{z - H^y 1 
to obtain 

= ~h L ( (z " Fp ° er ~ vr 1 ~ {z ~ H ° er)_1 ) dz = Qi + ° 2 + ° 3 

with 

Qi = -£- f(z- H^)- l V{z - H^-Hz, 

HIT Jcg 

Q2 = ^ J (z - HlJ- l V{z - H°J-W(z - H^)- l dz, 
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Notice that (1 - A)Q 3 £ 6 2 by (|5.12|) and the estimate 

H^-^J-'le. ^ ll^ 1 -^)" 1 !^ ll^)*||<C||F|| i6 

where we have used (|5.4p and Lemma [3] 

Let us now prove that (1 — A)Qi <E ©2- First we notice that 

[ (z- H^y^V^iz - H^dz 

= I (z-H a pc r 1 (l O pC r) ± V( 1 Q pc ^(z-H^ 1 .)- 1 dz = 
by the residuum formula. Then we have 

(1-A) / {z-HlJ-^V^l^iz-Hl^dz 

which belongs to ©2 by Lemmas [3] and [51 The proof is the same for Q^. □ 

5.5. Proof of Theorem [2| existence of a minimizer with chemical poten- 
tial. Let (Qn)neN be a minimizing sequence for (|3.1ip . It follows from (|3.10[) that 
(Qn)neN is bounded in Q. By Proposition [TJ (pQ n ) n eN is bounded in C n L 2 (R 3 ). 
Up to extraction, we can assume that there exists Q in the convex set /C such that 

i) Qn Q and |V|Q n |V|Q weakly in 6 2 ; 

ii) l< r - ^I 1/2 Q+ + I< r - ^l 1 / 2 - |< r - e^|V2Q++ [jff er _ £f| i/2 5 
|ff p ° or - e^Q--^ _ £F |V2 ^ _ e^Q-,^ _ ep \i/2 
for the weak-* topology of 61; 

iii) PQ n -^PQ inCnL 2 (K 3 ). 

Recall that Si is the dual of the space of compact operators [HI Thm VI. 26]. Thus 
here A n — 1 A for the weak-* topology of 61 means Ti(A n K) — ^ Tr(Aff) for any 
compact operator K. 

Then, as •) defines a scalar product on C, 

D(pq -v,pq~v) < liminf D(p Qn - v,p Qn - v). 

Now since Q+ + > 0, |if por — e F | 1 / 2 Q+ + |-ff pcr — ^f\ x ^ 2 is also a nonnegative operator 
for any n. Thus Fatou's Lemma [3D] yields 

Tr(|H p ° cr -e F | 1/2 Q ++ Kr -£f| 1/2 ) < liminf Tr(|# p ° er - e F \^ 2 Q++\H pcl .- e F \^). 

The same argument for the term involving —Q > yields 

£ V (Q) - e F Tr (Q) < liminf (f(Q n ) - e F Tr (Q„)) = , 

n — >oo 

i.e. Q £ K, is a minimizer. 

The proof that Q satisfies the self-consistent equation (|3 . 1 2[) is classical: writing 
that £ v ((l - t)Q + tQ) > £"{Q) for any Q £ K x and t £ [0, 1], one deduces that Q 
minimizes the following linear functional 
(5.13) 

Q£K,^ F(Q) := Tr(|tf p ° cr - e F \^ 2 (Q ++ - Q~~)\H° eT - erf' 2 ) + D(p Q - u,p Q ). 
Notice that when Q £ JC r C JC, one has 

F{Q) = Tio((Hq - e F )Q) 



LOCAL DEFECTS IN PERIODIC CRYSTALS 



21 



where we have used the definition of pq in Proposition [T] to infer 

D(p Q - v,p Q ) = Tr (((p Q -v)*\- r 1 ) g) , 

since (pq — v) -k \ ■ | _1 € C when pq — v E C. Minimizers of the functional (|5.13[) 
are easily proved to be of the form ()3.12p . They belong to /C r by Proposition [H □ 

5.6. Proof of Proposition [3j the value of E°(q). Clearly 

(5.14) E°(q) > inf{Tr (ff p cr Q), Q G K, Tr (Q) = q} := E°(q) 

since D(pq, pq) > 0. It can be easily proved that 

when q > 
E^g when q < 0, 

see, e.g., the proof of Lemma 13 in [T2"] . 

Thus it remains to prove that E°(q) < E°(q) which we do by a kind of scaling 
argument. Let us deal with the case q > 0, the other case being similar. We can 
assume that ^z+i ~ mm £er* Az+i(£) = Az+i(£o) since each A n (£) is known to be 
continuous on T* . For simplicity, we also assume that £o is in the interior of T* (the 
proof can be easily adapted if this is not the case). Let us denote by uz+i (£, •) € £| 
an eigenvector associated with the eigenvalue Xz+i{£.) for any £ E T* . It will be 
convenient to extend it on R 3 x R 3 by uz+i(S,,x) = when £ E R 3 \ T*. Since 
H° er u z +i(S,,x) = \z+i{i)u z+ i{i,x) for any £ E T* , it is clear that 

sup \\Au z +i(£,-)\\ L 2 (r) < oo, 

i.e. uz+i G L°°(R 3 ,i/ 1 2 oc (]R3^ and g £°°(R 3 x K 3) by the z 3 -periodicity 

(resp. the 27rZ 3 -periodicity) of e~^' x uz+i(£,,x) with respect to x (resp. to £). 

Consider now a fixed space V of dimension d — [q] + 1, consisting of functions 
/ : R 3 — * C with support in the unit ball B(0, 1) of R 3 . For any A > 1, we introduce 
the following subspacc of L 2 (R 3 ): 

W X := \gx ■= A 3 / 2 / /(A(£ - &))«z+iK, ■)<£> 



which has the same dimension as V by the properties of the Bloch decomposition, 
when A is large enough such that the ball S(£ojA _1 ) is contained in T* . Noting 
that for any gx E W\ arising from some f E V 



\gx(x)\ < A- 3 / 2 / |/(e>z+i(£o + \-^',x)\d? 

< \- 3/2 \\uz+l\\ L - {R3xm f \f(ZM, 

Jb(o.i) 

L»(E 3 xE 3 ) II/IIl 2 (R 3 ) ' 



3(0,1 

< CA- 3 / 2 | Uz+ i|^ fR 3 xRai |/| 
we deduce by interpolation that 



W2<p<oo, lim sup ||s| iP(R3) = 0. 

A^OO geWx 

ISIU2(R3)=1 

By construction one also has for any fixed f E V with associated gx E W\ 



* A — >oc 



0. 



(g x , (ff p ° or - ^ z+1 )gx) = A 3 J ^ |/(A(£ - £ ))| 2 (A z+ i(0 - £ Z+1 R 

Take now an orthonormal basis (<fii, ip^ +1 ) of Wx and introduce the operator 

[?] 

o a :=X;i^x^i + (?-[g])i^x^i. 



n=l 
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By construction 7p er <£ A , = f° r an y n = 1, [q] + 1 and Tr (Q A ) = Tr(Q A ) = g, 
thus Q x £ JC satisfies the charge constraint. Then 

[q] i v 

Tr (!/ p cl .Q A ) = £ (^,# p Vn) + (<Z M)(^ ]+1 ,if p ° er ^ ]+1 ) - £ z+1 <7 

n=l 

as A —>■ oo. Besides, (pQ\)\>i is a bounded family in i? 1 (K 3 ) which converges to 
in L P (R 3 ) for any p > 2. By the Hardy-Littlewood-Sobolev inequality [20], one has 

ii ii 2 

(5-15) D(p Q x,p Q x) < C\\pqx\\ l6/5 

which implies D(pgx, pgx) — > as A — > oo. Eventually £°(Q X ) — > £ (g) and 
Proposition [3] is proved. □ 

5.7. Density of finite rank operators in JC. This section is devoted to the 
generalization of results in [121 Appendix B] concerning the density of finite rank 
operators, which will be useful for proving Theorem [3] 

Lemma 8. For any Q £ JC there exists an orthogonal projector P such that P — 
7per ^ an d a trace- class operator S £ Q such that < S < 1, [P, 6] = and 

Q = P- Tper + S. 

Proof. This is an easy adaptation of O Lemma 19]. □ 

We denote for simplicity Sj+ = (1 - 7 p , cr )L 2 (K 3 ) and = 7 p ) cr L 2 (IR 3 ). 

Proposition 4. Lei P oe an orthogonal projector on L 2 (M. 3 ) such that Q = P — 
7p C1 . £ /C. Denote by (/i, /at) € (#+ n H 1 (M. 3 )) N an orthonormal basis of Ex = 
ker(P- 7 er -l)=ker( 7 er )nker(l-P) and by (gi,...,g M ) £ («- Hif 1 ^ 3 )^ an 
orthonormal basis of E_\ = ker(P — 7 p er + 1) = ker(l — 7 p Cr ) H ker(P). TTien i/tere 
eaiisi an orthonormal basis («i)i>i C fj+niP(R 3 ) of {E{) inSj+, an orthonormal 
basis (ui)i>i C £j_ fl P^IR 3 ) o/ (P-i) 1 - in £)_. and a sequence (Aj)j>i € 
swe/i £/ia£ 

(5-16) p = 2^|/„)(/„| + 2^ j— ^ . 

n— 1 2—1 1 

/, 17 x -, p f i w i , f K - KUj){Vj - \Uj\ 

(5-17) 1-P= 2^\9m){9m\+2^ ' 

m—1 i—1 1 

Additionally J2i>i (II Vm *IIl 2 + ll^illls) < °°- 

Proof. This is an obvious corollary of [TH Theorem 7] . □ 

Corollary 3. Let Q £ JC. Then there exists a sequence {Qk}k>i of finite rank 
operators belonging to /C r such that \Qk — Q\\q ~~ * as k — > oo and /or any k > 1, 
Tro(Q fc ) =Tr (Q). 

Proof. Taking A,; = for i > fc in the decomposition of Proposition |H one can 
approximate P by another projector P^ such that P& — 7 p G1 . — > P — 7 p er in Q as 
fc — > oo and P& — 7 p er is finite rank: 

N M k , 2 

(5.18) P k - 7 ° er = Y, l/n)</n| ~ E +E IT^O^I " !«*><"*!) 

n— 1 m—1 -i—1 * 

A 
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It then suffices to approximate each function in (|5.I8[) by a smoother one, for 
instance by defining for e -C I, iii ■= \\R e Ui\\ L2 R e Ui and orthonormalizing these 
new functions, where R e was defined previously in Equation (|5.8[) . 

Then for any Q = P — 7p Cr + 5 of the form given by Lemma [HI it remains to 
approximate 5 by a finite rank operator 8k such that [Pk, 8k] — 0, which is done in 
the same way. □ 

5.8. Proof of Theorem [3| existence of minimizers under a charge con- 
straint. The proof of Theorem [3] follows ideas of [12]. The proof that any mini- 
mizer solves (|3.12p is the same as before and will be omitted. 

Step 1: Large HVZ-type inequalities. Let us start by the following result, 
which indeed shows that (b)=>(c): 

Lemma 9 (Large HVZ-type inequalities). Let Z G N\ {0}, v G C and assume that 
(Al) holds. Then, for every q, q' G R, one has 

E»{q)<E»{q-q') + E Q { q '). 

If moreover there is a q' ^ such that E v (q) = E u {q — q') + E°{q'), then there is 
a minimizing sequence of E u (q) which is not precompact. 

Proof. Thanks to Corollary [3l the proof is exactly the same as [TZJ Prop. 6]. □ 

Step 2: A necessary and sufficient condition for compactness. The fol- 
lowing Proposition is the analogue of [T2l Lemma 8]: 

Proposition 5 (Conservation of charge implies compactness). Let Z G N \ {0} ; 
v G C, q S R and assume that (Al) holds. Assume that {Q n )n>i is a minimizing 
sequence in K x for p. lip such that Q n — Q G K for the weak-* topology of Q. 
Then Q n — > Q for the strong topology of Q if and only i/Tro(Q) = q. 

Proof. Let (Qn) r i>i ^ IC r be as stated and assume that Tro(Q) = q. We know from 
the proof of Theorem [2] that 

(5.19) £ V {Q) < lim £ v (Q n ) = E v (q), 

n — >oo 

hence Q G K. is a minimizer of E v {q). Therefore Q satisfies the equation 

Q = X(~eo,t F )(H Q ) - 7pcr + s 

for some tj? e (S^, ^z+i) anc ^ where S is a finite rank operator satisfying < 5 < 1 
and Ran(5) C ker(i?Q — e^). In particular Q G /C r by Proposition [2] We now 
introduce 

P ■= X{-oo,e F )(H Q ), P' := X(e P ,oc)(H Q ) and tt := X{€ F }(Hq). 
Let us write 

£ v {Q n ) = £ V {Q) + Tr (H Q (Q n - Q)) + l -D{p Qn - p Q ,p Qn - p Q ). 

Now using [TUJ, Lemma 1] and the hypothesis Ti'o(Q„) = Tro(Q), we obtain 

Tr (H Q (Q n - Q)) = Tr ((H Q -e F )(Q n -Q)) = Tvp((H Q -e F )(Q n -Q)) 
= Tr(\H Q - e F \(P'(Q n - Q)P' - P(Q n - Q)P)), 

where we recall that by definition Ti P (A) = Tr {PAP + (1 - P)A{\ - P)). We have 
-P - 6 <Q n ~Q < 1 - P - 6 which yields 

(5.20) P{Q n - QfP + P'{Q n - QfP' < P'{Q n - Q)P' - P{Q n - Q)P 
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and in particular P'{Q n -Q)P' > and P(Q n -Q)P < 0, i.e. Tr Q (H Q (Q n -Q)) > 0. 
Since we know that lim^^co £ v (Q n ) = £ L '{Q) 1 we infer 

(5.21) lim Tr(\H Q - e F \P'(Q n - Q)P') = lim Tr(\H Q - e F \P(Q n - Q)P) = 

n — >oo n — >oc 

and from ([5~2"U)) 

lim TrQH Q - e F \P'(Q n - Q) 2 P') = lim Tr(\H Q - e F \P{Q n - Q) 2 P) = 0. 

n — >oo n — >oo 

On the one hand, it is easy to see that 

P\H Q - e F \P > cP\H Q - k\P and P'\H Q - e F \P' > cP'\H Q - k\P' 

for some small enough constant c > and some k ^ (j{Hq) close enough to e F . On 
the other hand, the weak convergence of (Q n ) and the fact that tt is a "smooth" 
finite rank operator imply that 

(5.22) lim Tr(\H Q - K\n(Q n - Q) 2 ir) = lim TtQHq - k\tt(Q„ - Q)n) = 0. 

n — >oc n — >oc 

It is then clear that this yields 

(5.23) lim Tr(\H Q - «|(Q„ - Q) 2 ) = 0. 

n — >oo 

As we have chosen k, ^ (t(Hq), we can mimic the proof of Lemma [3] and obtain 
that 

(5.24) Ci(l-A) < \Hq-k\ <c 2 (1-A). 

Hence (j5T23l) shows that (1 - A) 1 / 2 (Q„ - Q) -> in 6 2 . 
Writing now 

(5.25) (Q n - Q)— = P(Q n - Q)P + ( 7 ° cr - P)(Q n - Q) 7 ° er 

- (P - 7 p Cr )(Qn - Q)(P - 7per) + 7per(Q« " Q)(7per " P) 

and using (pT2"4"|) . ([5~2"T]) and ([5~2"2"|) , we easily see that - k| 1 / 2 (Q„ - Q) — \Hq — 
kI 1 / 2 -> and (1- A) 1 / 2 (g„-Q) — (1 - A) 1 / 2 ^ in 6x. The proof is the same 
for (Q n - Q)++. □ 

Step 3: Proof that (c)=>(b). We argue by contradiction. Let {Q n ) n >i Q K> be 
a minimizing sequence for E v (q) which is not precompact for the topology of Q. 
By the density of /C r in /C, we can further assume that each Q n £ /C r . The bound 
(|3.9P on the energy tells us that (Q n )n>i is bounded in Q. Then, up to extraction 
and by Proposition \5\ we can assume that Q n — Q € /C where Tro(Q) ^ g, and 
that pQ n pq weakly in C. We write Tro(Q) = q — q' with q' ^= 0. We now prove 
that 

(5.26) E»(q)>E»(q-q>)+E°(q>) 

which will contradict (c) . To this end, we argue like in the proof of [TSJ Thm. 3] : 
consider a smooth radial function \ with support in B(0, 1) such that < \ < 1 
and x = 1 m £?(0, 1/2); define \r( x ) '■= x( x /P)- Then let be i]r := ^/l — Xr. Let 
us introduce the following localization operators 

Y R := 7per»?H7^r + (7pcr) V^pcr^, ^ = ^ " **• 

Lemma 10. We have for all 3 < p < oo, 

(5.27) ||[^,7pV]|| ep =0(i?- 1+3/p ), ||Yb - r, R \\ 6p = 0(iT 1+3 /*>). 
Moreover \\X R -x R \\ G = O (ir 1+3 / p ). 
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Proof. By QOj and (£ 

which yields ||[%,7poJ|| ep < <? ||(1 - A)" 1 ^)^ < C \\Vm\\ L p h Y the Kato " 
Seiler-Simon inequality and following the proof of Lemma [5J Eventually we notice 
Yr — T)R = (l — 27p Cr ) [7p 0r ,?7jj] and thus (|5.27[) is proved since 7p er is bounded. 
The last inequality is a consequence of ||Yr[ < 1, IvrI < 1- ^ 



Lemma 11. One /ias 



(5.28) 



Yr, \H 



E |l/2 



Iff 



E |-l/2 



where we recall that E = (E^ + E z+1 )/2 is £/ie middle of the gap. 

Proof. We use the well-known integral representation of the square root [2] 

i r i //;;..,. -^i 



(5.29) 



,,n yil/2 _ / l-^pcr S l 

' pcr ' Wo S + |ff p ° er -S|Vi' 



Recall that Yr = 7p l or ?7fi7p 0r + (7per) X7 ?-R(7per) ■ For shortness, we only detail the 
estimation of the term involving (7p 0r )" L ?7fl(7pcr) ± - Using that \Hp Cr — E| commutes 
with (7PO!.)" 1 " and that (7p Cr )" L is bounded, we see that it suffices to estimate 



(5.30) ( 7p ° er ) ± 



Vr, 



|g p ° CT -E| 
s + |ff° er -s| 



\ H pei 



e|-v^ (7 ° y 



(7° 
V /pc 



2j S + |ff° er -E| 



[ m -A}\H° eT -V\-V 2 - 



\ /pc 



\H° - E 



■y/sds 



Then [tjk,-A] = (At/*) + 2(Vry fl ) • V, hence ||[ry R , -A] |ff° er - E|-Va| = Otftr 1 ) 
where we have used that V|-ffp er — E|~ 1//2 is a bounded operator by Lemma[TJ Then 
we use that | ffp er — E| > e for some e > to estimate the right hand side of (|5.30[) 
in the operator norm by 

[^,-AllIf^-Sr 172 x r 7 ^- = 0(JT 1 ). 

Jo l e + S J 



□ 



Notice now that YrQ u Yr £ K. for all i? > 1 (the same is true for XjiQ n Xn but 
we shall actually not need it). To see this, notice for instance that 



(5.31) (1 + M)(YrQuYr) = -(1 + |V|)|ff p ° cr - E^/ 2 ( [y fl , |ff p ° or - ei 1 / 2 

x |ff p ° or - El" 1 / 2 !^ - E^Q,^ _ y fl |^0 er _ ^Q n Y R 



which belongs to & 2 since |ff° er - E| 1 /2g n e & 2 an d [y H) |#0 cr - E| x / 2 ] \H° er - 
E^ 1 / 2 is bounded by LemmaHU The proof that (1 + | V|)(YrQ„Yr) ++ (1 + |V|) 
and (1 + |V|)(YrQ 71 Yr) __ (1 + |V|) are trace-class is similar. Eventually, the proof 
that — 7p er < YrQYr < 1 — 7p CI . is easy, using the equivalent condition (|3.4[) and 
the fact that 7p Cr commutes with Yr. 
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We are now able to prove (|5.26[) as announced. We write, following [TT] . 

(5.32) £"(Q n ) = Tr(X R \H° ei - £| 1/2 (Q++ - Qn~)\ H per - A 1/2 X R ) 

+ Tr(Y R \H° pcl . - n l '\Q+ + - Q--)\H° pcl . - n^Yn) + ETr (X R Q n X R ) 

+ T,Tr (Y R Q n Y R ) + ^D(p Qn - u, PQn - v) - ~D(v, v) 

where we have used that [7° er , X R ] = [j° er , Y R ] = to infer Tr (Q„) = Tr (X R Q n X R )+ 
T^o(YRQnY R ). Then, by Lemma [TT] and using the fact that (Q n )n>i is a bounded 
sequence in Q, we deduce that 

(5.33) Tr(Y R \H° pci - £|V*(Q++ - Q~~)\H° ei - Z\V*Y R ) 

> Tr(|tf p ° cr - S| 1 /2 Yr (Q++ _ Q--)Y R \H^ ei - SI 1 / 2 ) - C/R 
for some constant C > 0. Arguing similarly for the other terms, we obtain 
(5.34) 

£"(Qn) > £°(Y R Q n Y R ) + Tr( XR \H° er - S| 1/2 (Q++ - Q~~)\H° er - X^xr) 
+ ETr( Xfl (Q+ + + Q--)XR) + \o{p Qn - v, PQn -v)- \d{v, v) - C'/R 
for some constant C", where 

£°(Q) := Tr(|if p ° er - E|Va ( Q++ - Q~ )|# p ° er - E\V 2 ) + STr (Q). 
Recall (Proposition [3]) 

E°(q)=M{£°(Q), Qe/C, Tr (Q)=g}. 

Thus, using 

« = Tr (Q„) = Tr (Y R Q n Y R ) + Tr (X R Q n X R ), 
and the fact that q i— > E°(q) is Lipschitz by Proposition [3] (|5.34[) yields 

(5.35) £»{Q n ) > Tr(xR\H a pcl ~ S| 1/2 (Q++ - Qr)\K* - S| 1/2 Xi?) 

+ £Tr( Xii (Q+ + + Q"-)Xil) + - Tr( X fl(g++ + 

+ \D{PQ n ~ v, P Qn -v)- \d(v, v) - C'/R 
Let us now pass to the limit n — * oo. First we notice 

(5.36) liminf Tr(x R \H^ ei - S| 1 / 2 (Q++ - Q-"Wer " S| 1/2 x fl .) 

> Tr( Xii |J? p ° er - E| 1 / 2 (Q++ - Q-)\H^ n 1/2 XR), 

(5.37) liminf D(p Q - ^, p Qii -v)> D(p Q ~v,p Q - v) 
by Fatou's Lemma and the weak convergence pQ n —^pq in C. Then 

lim Tr(x;?Q+ + Xfl) = ^(xrQ ++ Xr), lim Tr(xi?<3~~Xfl) = Tr(xi?Q Xr) 

n— »og n— s-oo 

which is obtained by writing for instance 

Tr(XRQt + XR) = Tr( X «(l + |V|)- X (1 + |V|)Q++(1 + |V|)(1 + M)" 1 **) 

and using that Xfl(l + | V| ) 1 is compact (it belongs to & p for p > 3 by the Kato- 
Seiler-Simon inequality) and that 

(1 + |V|)Q++(1 + |V|) - (1 + |V|)Q++(1 + |V|) 
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for the weak-* topology of &±. Thus, 
(5.38) 

E"{q) = liminf £»(Q n ) > Tr( XR \H° er - £| 1/2 (Q ++ - Q—)\H° er - ^\ R ) 

+ ETr( Xfl (g++ + Q—)xr) + E°(q- Tt( X r(Q ++ + Q~)xr)) 

+ X -D{ PQ -„, PQ -v)- ±D(v, u) - C'/R. 

Passing now to the limit as R — > oo, we obtain (|5 .26|) . This contradicts ([3]) and 
shows that (b)<=>(c) in Theorem [3J 

Step 4: Characterization of the q's such that (c) holds. Because q i— > E v (q) 
is a convex function, it is classical that the set / = {q e R, (c) holds} is a closed 
interval of R, see e.g. [21]. It is non empty since it contains Tro(Q) for any minimizer 
Q of E" p obtained in Theorem^ for any tp in the gap (S^, ^z+i)- Additionally, 
q i— ► E v (q) is linear on the connected components of R \ / and / is the largest 
interval on which this function is strictly convex. Let us now state and prove the 

Lemma 12. Let Z eN\ {0}, v G L^R 3 ) n L 2 (R 3 ), and assume that (Al) holds. 
Assume that Q\ and Q 2 are respectively two minimizers of E v (q\) and E v (q 2 ) with 
<Zi 7^ <?2- Then pq 1 ^ pq 2 and therefore 

E v {tq x + (1 - t)q 2 ) < £ v (tQx + (1 - t)Q 2 ) < tE v ( qi ) + (1 - t)E v {q 2 ). 

Proof. Assume by contradiction that pq 1 = pq 2 . It is classical that the operators 
Qi and Q 2 satisfy the self-consistent equations 

Ql =X(-oo,ei)(#Qi) -7per + 5 l) Qi = X(-oo,e 2 ) (Hq 2 ) - 7p 0r + 8 2 

where < 5k < 1 and Ran(Jfc) C kei(HQ k — e&) for k = 1,2. Necessarily ei and £2 
are in [Sj, otherwise Q\ and Q2 would not be compact, which is not possible 

since every operator of JC is compact. Since Hq 1 = Hq 2 has only a point spectrum 
in the gap, we deduce that if et € (S^, T>~^ +1 ), then necessarily 5k is finite rank. If 
e-k £ {^2, then it can easily be proved that at least 5k € ©i- Hence Qi and 

Q 2 differ by a trace-class operator: Q 2 = Qi + 5, Tr\S\ < 00. Now 7^ q 2 — q\ = 
Tr(5) = J ps which contradicts our assumption that ps — pq 1 — pq 2 — 0. The rest 
follows from the strict convexity of p t— > D(p, p). □ 

Corollary 4. There is no minimizer for E u (q) ifq^I, the interval on which (c) 
holds. Thus (a) implies (c) . 

Proof. Assume that there is a minimizer Q\ for some q\ £ I, for instance q± > 
max/ := q 2 . Applying Lemma fT2l to q\ and q 2 shows that E v (-) cannot be linear 
on [(72, qi] which contradicts the definition of I. □ 

5.9. Proof of Theorem [4j thermodynamic limit of the supercell model 
for a perfect crystal. 

Step 1. Let us first prove that limsup yS-^sc l — ^per Starting from the Bloch 
decomposition 

1 

(2V) 3 

of 7p er it is possible to construct a convenient test function j SCy L for (|4.2| as follows: 

y > j e ^ Z 3 nr . we+i — 1-, — E — ) / 



7pcr - 775^ 
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with T) = if L is even and 77 = 1/2 if L is odd. It is indeed easy to check that 7 sc ,l 
is in V sc ,l and satisfies Py ao L — P-y° eT - In particular, 



Ja l 

and, since both p 7 o^ and /i pcr are Z 3 -periodic, 

D G L (Py sc . L - Mper, Py sc , L ~ Mpcr) = L 3 D Gl [p^ - ^ per , p 7 o - /X pcr 



Besides, 

^^s-caj (-\ A ^) = (2^)3 / r . (4 A W-'s) * 

"(2^ S L „ (e-e')-^cr)Hv( 7 ° er ),)de'. 

It follows from the boundedness of J r » Tr L 2( r -,((l — A)(7p er )j) and from the 
inequality | — 2iV| < (1 — A) that the last two terms of the above expression go to 
zero, hence that 

Therefore lim L ^ +oc L" 3 £° c L {% c x) = £ pcr ( 7p cr) = ^pcr and consequently 
limsup ^/ S ° C . L < lim (j sc L ) = ^ 

Step 2. Let us now establish that liminf -^I^ c L > 1°. First, the existence of a 

L— >+oo ' P 

minimizer ji for (14. 2| and the uniqueness of the corresponding density p® c L follows 
from the same arguments as in the proof of [5J Thm 2.1]. Note that, by symmetry, 
p° c L is Z 3 -periodic. We now define the operator 7° c L as 

1 \ ~» * 

7sc,L = 2.^ T klLTk- 

By simple periodicity arguments, it is clear that t^lth is also a minimizer for 1)4. 2p 
for all fc e Z 3 . By convexity, so is 7g C L . Besides, 7° c L commutes with the transla- 
tions Tfc for all k G Z 3 so that its kernel 7° c L (x, y) satisfies 

V(x, y, z) G M 3 x M 3 x Z 3 , + z, y + z) = 7 ° c , i 0r ! 2/)- 

The optimality conditions imply that j®c l can be expanded as follows 

£e^z 3 nr» fc>i 

where for any £ G ^Z 3 n T*, (w^)fe>i is a Hilbert basis of £p er (r) consisting of 
eigenfunctions of the self-adjoint operator on Lp Gr (r) defined by 

~A - • V + (p° c , L - Mper ) * r G x + 
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associated with eigenvalues Af (£) < (£)<■• • The occupation numbers n£ ^ are 
in the range [0, 1] and such that 

r* E E<« = z - 

^ e ^LZ3 nr » k>l 

Lastly, there exists a Fermi level £^£l such that n^c = 1 whenever (£) < e^. 
and = whenever A£(£) > e^.. One has 



(5-39) p0 cL = ^_ g cL{j cL) ^^_ £ £l|i|| HV + ^y^ (r) 

£e^z 3 nr* fc>i 
1 

2~L 3 

We now introduce 



+ ^^(/M " Mpcr,P 7 o ~ Mpcr). 



(2tt) 3 

where /? L (f ) is the unique element of ^Z 3 nr* such that £ e MO+h^ 2 ^"> ) 3 . 
It is easy to check that 7° c L G Pp er - Thus 7° c L € Vp er can be used as a test function 
for (El). Therefore 



(5-40) ^Al) > 7 per- 

As p^o = p^o is Z 3 -periodic, one has 

•sc,L 'sc,L 

(5.41) D Gl {py<> c L ~ Mpor, PjO c L - Mpcr) = J1 D Gl (P 7 ° c l ~ Mpor, Ay^ ^ - /i per ). 



Besides, 
(5.42) 



(2^ 1 *W (4 A W^) « - (2^ I g ^"(-^^^(dl^CD 



h E E^IIH V + ^ll^ r (r)+^ 



with 



^J^O (||(-iV + eKV K) ll! ?er( r) - IK-^ + fe(OK^(oll^(r)) 



r fc>l 



(27T) 

Putting (dHU) dnHU) together, we end up with ^l" c L + R L > /p er - As 



\Rl\< 



L 3 I L 2 



we finally obtain liminf -pr/ s ° c L > I^ eI - 

Step 3: Convergence of the density. A byproduct of Steps 1 and 2 is that 
(7sc i)ieN\{o} is a minimizing sequence for 7°. The convergence results on the 
density p~a = p^a immediately follow from the proof of 5, Thm 2.11. 
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Step 4: Convergence of the mean-field Hamiltonian and its spectrum. 

One has H® c L — i?p 0r = $l where $l solves the Poisson equation — A$£ = 
47r(/9 7 o ^ — p-fO ) on r with periodic boundary conditions. As it follows from Step 3 
that (p 7 o — p 7 o or ) converges to zero in Lp Cr (F), we obtain that <J>l converges to 
zero in Hp eT (T), hence in i°°(R 3 ). Consequently, 

(5-43) ||^ s c ,L-^p e r || < l*£li«-0 

as L — > oo. This clearly implies, via the min-max principle, that 
sup sup |a£(0 - A n (0| < — » o 

n>l{6r* L^oc 

where (A£(£)) n >i,£ er * (resp. (A n (£))„>i, f er *) are the Bloch eigenvalues of iJ s ° c L 
(resp. if° er ). 

Step 5: Uniqueness of 7° c L /or large values of L. In the remainder of the 
proof, we assume that (Al) holds, i.e. that H® er has a gap. 

The spectrum of H® c L considered as an operator on Lp er (Ai) is given by 

™eN\{o}^e^z 3 nr* 

It follows from Step 4 that there exists some Lq £ N \ {0} such that for all L > Lq, 
the lowest ZL 3 eigenvalues of H® c L (including multiplicities) are 

U U tffc) 

\<n<Z £g^LZ3 nr . 

and there is a gap between the (ZL 3 )-th and the (ZL 3 + l)-st eigenvalues. As a 
consequence, 7° c L is uniquely defined: it is the spectral projector associated with 
the lowest ZL 3 eigenvalues of H® c L , considered as an operator on Lp Cr (Ai). 

Step 6. Let e F G (E^,E^ +1 ). For L large enough, 7° Cj£ = X(-oo,e F ] (H° C<L ) = 
X(~oo,o](Hsc.l~ £ f) as operators acting on Lp er (Ai). This means that Tsc.l satisfies 
the Euler-Lagrange equation associated with /° Cjijejr , see (|4.5|) . As the functional 
7 i— > £^(7) — eFTr L 2 o ( Ai )(7) is convex on V S c,Li 1sc,l ^ s a minimizer of this 
functional. Its uniqueness follows as usual from the uniqueness of the minimizing 
density and from the fact that is not in the spectrum of _ff s ° c L — ep. 

5.10. Proof of Theorem [5} thermodynamic limit of the supercell model 
for a crystal with local defects. We follow the method of [13]. As in the 

previous section, we denote by 7° c L the minimizer of l|4.2p , which is unique for 
L large enough and is also the unique minimizer of (|4.5[) . Let K-l be the set of 
operators Ql on Lp C1 .(A L ) such that 7sc,Z + Ql e 'Psc.l- In fact 

K L = {Ql € &x(Ll eI (A L )) | Ql = Q L , |V|Q L |V| € S^L^Az)), 

-7 s °c,l<Ql<1-7sc,l}. 
We introduce £^. tLt£F := S^c.l — £ FT^L 2 et (A L )- Let Q L e /Cl, one has 



' \ D Gl (pQl > PQl ) - efTr L 2 cr (At } (Q L ) - D Gh (u L , p^ L - p pcl ) + ^D Gl (y L , u L ) 
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Note that in the above expression, H® c L is considered as an operator on Lp er (Ai). 
Using Theorem HI this equality can be rewritten, for L large enough, as 

(5.44) 

Cc,l, £f (7 s °c,l + Ql) - £ Q sc ,l,,Mc.l) = Tr i?er(A£) (|tf s c>i - e F \(Q++> L - Q^)) 

+ \ D G L {PQ L ~ V LiPQl - V L ) - D Gl (v L , Pi o^ l - PpcI ) 

where we have set 

Ql + ' L = (1 - 7 s °c,l)Ql(1 ~ 7 s °c,l) and Q^ L = 7° sc<L Q L J° CtL . 
It follows from (15.44)1 that 

C,L, iF ~ 4°c,L, eF = inf {e^ l (Ql) - e F Tr L2paAAL) (Q L ), Q L £lC L ) 

- D Gl {v LiPi ^ l - Ppcl ) + \^D Gl (v l ,v l ). 

where 

E sc,l(Ql) - e F Tr L 2 er(Ai) (Q i ) := -D Gl ( PQli v l ) + \d Gl ( PQl , PQl ) 

+ Tr i?er(Ai) - eF |V»(g++.^ - Q, •'-,//■'., - e ,|V») . 

First, using ^ being in L 1 (M 3 ) n L 2 (M 3 ) and the convergence of $l = (p 7 o ^ — 
p 7 o ) *r Gi to zero in L°° (see Step 4 of the proof of Theorem H] in Section [B3f . 
we obtain 

If 1 

-D GL {v L ,po - Ppci ) + -D Gl (v l ,v l ) — ► - / v({ Pl o - Ppcl )*rG 1 ) + -D(v, v). 

I J R 3 P I 

Our goal is to prove that 

(5-45) lim £ = El 

where 

(5.46) E» FtL = inf {E^ L (Q L ) - e F Tr LlM (Q L ), Q L E /C L } . 

Step 1: Preliminaries. In the proof of (|5.45|) . we shall need several times to 
compare states living in L 2 er (Ai) with states living in L 2 (R 3 ). To this end, we 
introduce the map 

i L :L 2 (R 3 ) - L 2 C1 .(A L ) 

zGZ 3 

Notice that (il)* ■ L 2 ci (Al) — > L 2 (M. 3 ) is the operator which to any periodic 
function ip 6 Lp er (A/,) associates the function Ia^V G L 2 (K. 3 ). Remark that 
*l(*l)* = (Ai) whereas (zl)*«l = J-Ax,- Hence ix, defines an isometry from 

L 2 (A L ) to Lp er (A_[,). The equality = <£> is only true when 93 e L 2 (R 3 ) has 

its support in Al. When tp e i/ 1 (K 3 ) satisfies Supp(<y9) C Al, then one also has 
dxjL(cp) = ih{d Xi (p). 

Notice in addition that if A G 6i(L 2 cr (A L )), then {il)*Ai l £ &i(L 2 (A L )) C 
6i(i 2 (M 3 )) and 

Tr i?e r (Ai,)(^) = Tr L 2 (R 3 ) • 

Similarly if A S 6 p (L 2 cr (A L )), 

( 5 - 47 ) I^IspC^CA,)) = \( l L)*AlL\\ &p[L 2 m) . 
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Finally, we shall use that for any Z 3 -periodic bounded function /, i^f — fi^, 
where we use the same notation / to denote the multiplication operator by the 
function / acting either on L 2 oi .(Al) or on L 2 (R 3 ). Similarly, the operator H® c l — 
—A/2 + (p 7 o ^ — ^ pcr ) *r Gi can be seen as acting on Lp er (Ai) or on L 2 (R 3 ) and 

we use the same notation in the two cases. Then we have for any cp E Cq°(R 3 ) 
satisfying Supp(y>) C Al 

(5-48) K,LiLV = i L H° aCjL y. 

Notice that one can also define — iV on L 2 (R 3 ) or on L^ el {Ai) and we shall adopt 
the same notation for these two operators. We gather some useful limits in the 
following 

Lemma 13. Let be ip E L 2 (R 3 ) and tp E C£°(R 3 ). Then we have as L oo 

(1) {itTll, L iL^^lUi> m L 2 (R 3 ); 

(2) {i L yH^ Ll ^ L i L <p -> ff p ° er7p > m L 2 (R 3 ); 

(3) (i L )*A^ CiL i Lt p - A 7 or ^ m L 2 (R 3 ); 

(4) {i L )*{l + \V\)i L <p -»■ (1 + |V| V in L 2 (R 3 ); 

(5) (iL)*\H° CiL ~ e F \ l /H LV - |if° er - m L 2 (r3) /or flnj/ ^ £jr m 
the gap (£+,£ z+1 ). 

Proof. The operator (iL)*7s ci i_L being uniformly bounded with respect to L, it 
suffices to prove the first assertion for a dense subset of L 2 (R 3 ) like C^°(K. 3 ). Hence 
we may assume that ip = tp E Cg" (R 3 ) . 

Let K be a compact set in the resolvent set of H® CI . We are going to prove that 

(5.49) lim (i L )*(z - H^ L )-H L <p (z - ff° er )"V 

in L 2 (R 3 ), uniformly for z G AT. To this end, we first notice that by Theorem^] K 
is contained in the resolvent set of H® c L for L large enough and thus 

\\{ z ~ -^sc,l) ~ ( z ~ ^per) IIb(L 2 , (A l )) — ^("^0 ||-^pcr — -^sc,i||g(£2 (Al)) > ^ 

by (|5.43|) . Hence it suffices to prove (|5.49p with _ff " c L replaced by H pel (seen as an 
operator acting on L 2 6I .(Al)). Then we use its Bloch decomposition (detailed in 
Appendix, see (|A.2|) ) and compute, assuming L large enough for Supp(</?) C Al, 

{i L )*(z-H° et )- 1 i L ip(x) = V V — L ( [ e n (k,-)ip)lA L (x)e n (k,x), 

2 L~ 3 f 

It is then easy to see that {il)*{z — H^^^Hup — ^ (z — H® )~V weakly in L 2 (R 3 ) 
(one can take the scalar product against a function tp E Qj°(R 3 ) to identify the 
limit) and that ||(i L )*(z - ijO er )~ 1 i L ^|| i2(R3) -> ||(z - H^)- 1 ^^^, yielding 

the strong convergence in L 2 (R 3 ). 

For the proof of (1), it then suffices to choose a curve around the first Z bands 
of Hp Cr and use the above convergence of the resolvent in the Cauchy formula. 
Assertion (2) is an easy consequence of (1) and (|5.48[) . Indeed, by Theorem [H we 
know that lim^oo \\H ( J c L - Hper\\ B ( L2 ( R3 ^ = °- Since ( i i)*7 s ° c ,L i i is bounded, this 
implies that for L large enough such that Supp(</?) C A^ 

\\(iLT7° sc , L (H° C}L i L - lL H° p J<p\\ L2(m = \\(il)*1 ( : c , l il(HI l - H° p J<p\\ L2(R3) - 
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where we have used (|5.48p . Then we notice that H® cr tp E L 2 (R 3 ). Hence (1) implies 
that limL^oo || ((zl)*7 s ° c ,z,*£ ~ lper) H peM\ L ^(V) = °- The argument is exactly the 
same for the third assertion (3) . Assertion (5) can be proved in the same way, using 
(|5.48p and the integral representation of the square root ()5.29p . 

Finally, it remains to prove that (4) is true, which is done by computing explicitly, 
for L large enough such that Supp((p) C Al, 

(i L )*(i + M)i L v= ^^(i + \k\mky k - x tA L (x), 
\{i L y{\ + |v|)i^ii 2(R3) = J2 ^l(i+|fc|)^)l 2 -l^oo |(i + |v|)^ii 2(R3) . 

The strong convergence is obtained as above. □ 
Lemma 14. Let V <E C£°(R 3 ). We have as L -> oo 

(5.50) (t£)*(l - A)- 1 z L (^ L -» (1 - A)" 1 ^ 

(5.51) + |v|)- 1 i i (y)(i + M)-H L - (i + |v|)-^(i + M)- 1 

strongly in & 2 {L 2 (R 3 )). 
Proof. For L large enough, we have 

\\(i L y(i - ArH L (v) lL \\ 62{L2{m) = ||(i - Arw(v)\\ 62{L2pM) 

< \\vh 1 

- (2tt) 3 / 2 

\ Alt Mj ~ 

which shows that (?l)*(1 — A)~ 1 i[ j (V)iL is bounded in 6 2 (£ 2 ( 

(5.52) to £ -^§y=l * ,(p) = d + b| 2 )- 1 . 

Arguing as in the proof of the fourth assertion of Lemma 1131 we can prove that 
(|5.50p holds in the strong sense, hence the convergence holds weakly in &2(L 2 (M. 3 )) 
towards (1 - A^V . Now 

T II/- \*M A\-l- nr\- II lrlL 2 (R3) |5Il2( R 3) ,, A S-l T r|| 

^lim \\(i L ) (1-A) tL(V)i L \\ e2(L2{u3)) = ( 27r ) 3/2 ; = ^|| e2(i2(R3)) 

and the limit holds strongly in & 2 (L 2 (R 3 )) . 

The argument is the same for (|5.51|) . noticing that 

WlT(i + |v|)-V(y)(i + |v|)-V|| 2(i2(R3)) 

= Tr i2 (Al) ((1 + |V|)- 2 iL (V)(l + M)~ 2 i L (V)) 




(2n)- 3 / 2 ff \h L (x~y)\ 2 V(x)V(y)dxdy 

(i + M)-Mi + |v|)-f e2(i2(H3)) 



* L — ► oc 

where we have used that 

fo~\3/2 

Ak 



^ (2tt) 3 / 2 . 



converges to the Fourier inverse T 1 (h) of h(p) = (1 + \p\) 2 , strongly in L 2 oc (R 3 ). 

□ 
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Step 2: Upper bound. We prove here that limsup^^^ L < E" F . Let e > 0. 
Using Lemma [SI Proposition [H Corollary [3J and the notation therein, one can find 
a finite rank operator Q € /C r such that 

(5.53) E" p < £ V {Q) - e F Tr (Q) < E» p + e, 
of the form 

-1 -1 k A2 

(5.54) Q= E \ v ™){v m \- E KK U "1 + E i + V? (\ v i)( v i\ - K)("*l) 

m=-M n=-N i=0 1 

+ E 1 i + \ v i)( u i\) + 5 ' with ^^/J^'hjX^jl- 

*=o 1 + A * j = l 

Let < r\ « 1. It is possible to choose a family of orthonormal functions u 7 ^, v^, 
such that 



(5.55) \\vP n - u n \\ H 2 < r), \\v^ - v m \\ S 2 < T], \\w] - Wj\\ H 2 < 77. 

for all n = —N...k, m — —M...k and j = 1...J. Let us define the Gram matrices 

which, by (I5.55|) satisfy = Id M +k+i + o(l)^^ and Si = Id N+k+1 + o(l)^^ - 
We also introduce the orthogonal projector IT' on Span{7p Cr u^, (1 — 7p er )u^} and 
define (S*)ij := ((1 - W)w^v^). Clearly 52, = /dj + o(l)^ . 

nnrmnl cv«tflm in 7^ | 



Now we introduce a new orthonormal system in Lp Gr (Ai) 



: = E ( S '-,L 2 )i,n7 S c 1 i*i<) : = E (' S +, 1 / 2 )i,m(l - 7 s ° c ,L)^m 



=-JV m=-Af 



(S-,L)i,j = (iscXlLul ,i L uJ) L 2 er(Ai) , = ((1 -7s cx)*L^^i U j) i 2 (!r(Al/) - 

Notice that by the first assertion of Lemma I13[ hin^^oo <S±.l = S±. Finally, we 
introduce the projector 11^ on Span(u^ i; i/^ ^) and define 

J 

W ],L := X]^.i 2 )^^ ~ n i)*i w ^ ( S w,L)i,j = ((1 - ^L)iLwV,i L W^}. 

e=i 

We now define a state in JCl by 

Ql= E KU)«J- E I<l)(<J+E t ^(Kl)(<J-I<l)KJ) 



=-M n=—N i=0 

+ E TTA2 (Kl>KlI + Kl)«J) + X>Kl>KJ- 

i=0 2 j=l 

By Lemma fT3l we have 

(5.56) (ii)*< i£ -*l^oo <, («l)*<L,l ^l^oo % and {i L )*w^ L ^l^oo 
in L 2 (K 3 ) n L°°(R 3 ), where the limits are defined by 
k k 

V'i ■ ^ ] i^—)i,n 7por u ni ■ ^ ] (^+)i,m, (^ — 7pcr) u mJ 

n——N m=— M 

J 

.; 

£=1 
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By Lemma [TUl we know that for any fixed ip, ip € C^°(R 3 ), 

iil^o ~ e F)lsc.L i Lf, (ii)*7sc,i^^) = ((^per ~ ^hpexV, 7per^) ■ 

Hence, inserting the definition of Q\ in the kinetic energy and using the convergence 
of the Gram matrices, we obtain 

L lim Tr^ er(Ai) ((H^ L - e F )Ql) = Tr L2(R3) ((ff p ° er - e F )& 
where Q v is defined similarly as Q\ but with the functions (uJJ, vj^, instead of 
Let us now prove that 

L^yo Dgl ^ ,P ^ = D (PQ" ' p Q^- 

The convergence (|5.56j) implies that 1a l Pq v converges to Pq v in particular in 
I^QR 3 ) n L 2 (R 3 ). Notice the definition of Dq l (•, •) implies that 

(5.57) Vp e i^ cr (A L ) n L per (A L ), D Gl (p, p) < C (||yo||^ er(Ai) + \\p\\l leA A L )) 

for a constant C independent of L. Let us now write pq^ = p%,L + Pi,l where 
pi y L is the periodic function which equals ^b{q,l/4)Pq v l on At- The convergence of 
1 Al p Q i L towards Pq v in L^M 3 ) n L 2 (M. 3 ) and ([537| give that 

lim ||p2,l|| l i (Ar) = Hm |p 2 ,i|| i2 fArl = r lim d g l (P2,l, Pi,l) =0. 

L^OO perV ^ L— >00 PcrV ^ V L— >00 

Hence, it remains to show that 

(5-58) lim D GL (p hL ,p 1:L ) = D(pq v ,P£ v ). 



To this end we use the estimate [32] 

f 



sup 



G L (x) 



to obtain 



Dg l {pi,l,Pi,l) = // G L (x-y)pi, L (x)pi, L (y)dxdy = D(l Ai/ 9i jL , 1 Al Pi,l)+0(L x ) 
J J (At) 2 

where we have used that ^ Al ^ is uniformly bounded and that x — y G 

for any x,y £ _B(0, L/4), the support of pi t L- The convergence of 1a l Pi,l towards 
Pq v in L 1 (M 3 ) n L 2 (R 3 ) then proves (|5.58|) . Using the same argument for the term 
D G L (Pol , v h) we obtain 



lim ( -D Gl (pq*,,v l 

L — >oo 



Finally 



(5.59) lim E^ L (Ql) - e F Tr L?or(Ai) (Q£) = - e F Tr(Q"). 



L — >oc 



Passing to the limit as rj — > using (|5 . 55|) and the convergence of the Gram matrices 
S± and S^, we eventually obtain 

limsup^L < £"{Q) - £FTr(Q) < E» p + e. 

L — >oo 
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Step 3. Lower bound. We end the proof by showing that lim inf z,_,oo E" L > 
E^ F . As Ai is bounded for any fixed L, the existence of a minimizcr Ql of 

inf {E^ L (Q L ) - e F Tr i?cr(Ai) (Q i ), Q L e K L ) 

is straightforward. In addition, the spectrum of H® c L , considered as an operator 
on Lp er (Ai), being purely discrete and bounded below, Ql is finite rank. 

Using (14. 4p and reasoning as in the proof of Lemma [1] (see Section I5.1.1j) , we 
prove that there exists a constant c > (independent of L) such that \H® C l ~ € f\ > 
c(l — A) on ip Cr (Ai), for L large enough. The following uniform bounds follow 
from Step 2: 

(5.60) Tr L , cr(Ai) (|ie L - e F \ 1/2 {Q + L + - L QZ^Kl - *f| 1/2 ) < C, 

(5.61) Tr LlM ((l + |V|)(Q++- L - Qr' L )(l + |V|)) < C, 
(5-62) Tr i?et(A£) ((l + |V|)Q|(1 + |V|)) < C, 
(5-63) D Gl ( PQl - v l , PQl - v L ) < C, 

with C independent of L. 

Consider now the sequence of operators Ql '■— (ihYQ^L acting on L 2 (R 3 ). It 
is bounded in & 2 (L 2 (R 3 )) by (j5"^2"j) and since Tr L 2 (R3) (Q 2 L ) = Ti L2 ex (A L ){Ql) h Y 
(|5.47|) . Hence Ql weakly converges, up to extraction, to some Q £ &2(L 2 (R 3 )). 
Similarly, the Hilbert-Schmidt operator Rl ■= (iL)*Qi(l + |V|)«L weakly converges 
up to extraction to some R in & 2 (L 2 (M 3 )) . Let <p and if) be in C^°(IR 3 ) and assume 
that Supp(</?) U Supp(-0) C A L . Then 



<(u)*ol(i + ivi)^,^> i3(R3) = (QUii.ni + w)iL<p,r /V!ii 

-^L^oo (Q(l + |V|)cp, ^) jL 2 (R3) , 

where we have used that Ql Q weakly in 62 and that (%l)*(1 + |V|)ir,¥> — > (1 + 
\V\)<p strongly in L 2 (M. 3 ) by the third assertion of Lemma [T3l Hence Q(l + |V|) = 

Re 6 2 (L 2 (K 3 )). 

Similarly, define the operator Sl '■— (iL)*Q^' L ^L which is nonpositive and 
yields a bounded sequence in 6i(L 2 (IR 3 )) by (|5.6ip . Up to extraction, we may 
assume that (Sl) converges for the weak-* topology to some S S 6i(L 2 (IR 3 )). To 
identify the limit S, we compute as above for tp, ip G L 2 (IR 3 ), 

(SlV,iP) L 2 {R3) = <(*i)*7 s °c,L ( 9L7 s c,L*iV J V') L 2 (R 3 ) 

Using now the first assertion of Lemma [TBI we obtain limi^oo (Sl<P,iI>) l 2 (k 3 ) = 
(Q tpyi/j) i2(R3)- Hence Q = S S Si. The same arguments allow to conclude 
that in fact, Q € /C. 

Now, let T L := (i£)*l#sc,i - zf\ 1,2 Q~l~' L \ H s^l - e F \ 1/2 i L which also defines a 
bounded sequence in 6i(L 2 (E 3 )). Up to extraction, we may assume that Tl T 
for the weak-* topology of ©1. Arguing as above and using Lemma fl3l we deduce 
that T = \H° eI - e F | 1/2 Q" \H° cr - e F | 1/2 . Now, Fatou's Lemma yields 

lim inf Tr i?K(Ai) (\H^ L - ^(-Qr'^W^L - ^| 1/2 

= lhninf Tr i2(R3) (-T L ) > Tr i2(R3) (\H° pm . - e F \^ 2 (-Q—)\H° pcl . - e F \ 1 '^ 



LOCAL DEFECTS IN PERIODIC CRYSTALS 



37 



This proves that 

liminf Ti L 2 {Al) ((H° c<l - e F )Q L ) > Tr (i/ p ° or Q) - e F Ti (Q). 

We now study the term involving the density pq l . First, following the proof of 
Proposition [T] and using the bounds (I5.60|) - (|5.62j) . we can prove that there exists 
a constant C such that for all L large enough — ^- Hence, up to 

extraction, we have ^a l Pq l — 1 P weakly in L 2 (R 3 ) for some function p 6 L 2 (R 3 ). 
We now introduce an auxiliary function pl E L 2 (R 3 ) defined in Fourier space as 
follows: 

where for any k S (27r/L)Z 3 \ {0}, B k ■= B (k + 10 ^ k \ , tjjt;) which is chosen to 
ensure that l/\k'\ < l/\k\ for any k' S B k , and B := B (0, j^). 
Notice that pl is bounded in L 2 (R 3 ) as we have by definition 



/ p\=\ Xpl?= £ Mpq,)\ 2 = I 

JM 3 JR 3 J A 



ke^-z 3 



PQl 

Al 



On the other hand (up to extraction) p^ — 1 p weakly in L 2 (R 3 ), the same weak 
limit as t\ L pQ L . This is easily seen by considering a scalar product against a fixed 
function ip £ C^°(R 3 ). Now by the choice of the balls B k , we also have for L 3> 1 

D(p L , PL ) = 4rr I ^^ dki < D Gl { PQliPQl ) < C. 

Hence, up to extraction we may assume that pz, — p weakly in C. Using the 
regularity of 9, we also deduce that 

(^D Gl ( PQl ,v l ) + ^D Gl (pq l ,pq l ) \ >-D(p,v) + ~D(p,p). 

What remains to be proved is that p = pq where Q is the weak limit of (?l)*Ql*l 
obtained above. This will clearly show 

lraiinf K rt L > £"(Q) ~ CFTro(Q) > E» F 

and end the proof of Theorem We identify the limit of 1a l Pq l using its weak 
convergence to p in L 2 (R 3 ). 

We start with Pq++,l and write, fixing some V € C^°(R 3 ) and assuming L large 
enough for Supp(V) C A^, 

/ Pq + l + - lV = Tr Ll er (A L )(Qi + ' L iL(V)) = Tr L 2 (R 3)(A L B L ) with 

A L := ( lL )*(l + |V|)Q+ + ' i (l + |V|) iL , B L := (u)*(l + |V|)- 1 i i (F)(l + |V|)- 1 i i . 

The sequence (A L ) is bounded in ©i(L 2 (R 3 )), hence in & 2 (L 2 {R 3 )), by (1531) and 
converges (up to extraction) towards (1 + |V|)Q++(1 + |V|) weakly in 6 2 (£ 2 (R 3 )) 
(we proceed as above to identify the weak limit using the fourth assertion of 
Lemma [T5|) . By Lemma Bl converges towards (f + | V|) _1 y(l + |V|) _1 strongly 
in 6 2 (L 2 (R 3 )). We thus obtain 

lim / P q ++.lV = Tr L 2 (R 3)(<5 ++ U) = / p Q ++V. 
l^ooJ Kl *l J r3 

Likewise, it can be proved that the weak limit of Pq--,l is Pq--- 
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Let us now treat Pq+-,l (the other case Pq-+,l being similar). Following the 
proof of Proposition [U we write 

/ PQi-*V = Tr L 2 or(Ai) (Q + l - L [i° sc , l MV)]) 

J Al 

= ~ dzTr LU{AL) (Q+- L (z H° c<L y l (Ai L (V) - i L (V)A)(z fl^)" 1 ) . 
We only detail the argument to pass to the limit in 
Tr^A.) (Qt' L (* ~ H^ L yH L {V)A{z ~ K.Lr 1 ) 

= Tr L?or(AL) (A(z - H^ L )^Q + L -' L {z H^ L )-\l A)(l - A)" 1 ^)) 

= Tr L2(R3) (C L (i L )*(l - A)-H L (V)i L ) 

with C L := (i L )*A(z - H% L )- l Ql~' L {z - H° cL )- x {l - A)i L . One has, up to 
extraction, C L A(z - H° CI )- 1 Q+-(z - H^y 1 ^ - A) weakly in 6 2 (£ 2 (R 3 )). 
To see this, one first remarks that Cl is bounded in ©2(£ 2 (R 3 )) and then identifies 
the weak limit by passing to the limit in (C^ip, for some fixed <p, i\> € Qf(R 3 ), 
using the uniform convergence of the resolvent for z G as shown in the proof 
of Lemma fT3l Then by Lemma IT4l we know that (zl)*(1 — A)~ 1 ii J {y)iL converges 
towards (1 — A) _1 V strongly in &2(L 2 (Al)), hence we can pass to the limit in the 
above expression, uniformly in z € . We conclude that 

lim / 1 AlPQl V = / p Q V 

for any V S C^°(R 3 ), thus p = pq. □ 

Appendix A. Proof of Theorem [TJ 

Our proof uses classical ideas for Hartree-Fock theories. See [321 Section 4] for a 
very similar setting. Let us consider a minimizer 7p er of 7p er (it is known to exist by 
5, Thm 2.1]). First we note that the periodic potential V peT := (p 7 o or — p pC r) *r Gi 
is in L 2 oc (R 3 ). Thus V peI defines a A-bounded operator on L 2 (R 3 ) with relative 
bound zero (see [13 Thm XIII. 96]) and therefore H° el = -A/2 + V per is self-adjoint 
on V(-A) = H 2 (R 3 ) with form domain H 1 (W 3 ). Besides, the spectrum of H° eT 
is purely absolutely continuous, composed of bands as stated in J3H Thm 1-2] and 
[271 Thm XIII. 100]. The Bloch eigenvalues \ k (£), k > 1, £ e T* are known to be 
real analytic in each fixed direction and cannot be constant with respect to the 
variable £. Hence the function 

C:^^IUer* | A fc (0 < «}[ 

k>l 

is continuous and nondecreasing on R. The operator H® cr being bounded from 
below, we have C = on (— oo, inf Afe(r*)) and it is known [MJ Lemma A-2] that 
lirriK^oo C(k) = oo. We can thus choose a Fermi level e F such that 

(A.l) Z = C{e F ) = Itf e T* I A fc (0 < e F }\ . 

k>l 

Considering a variation (1 — i)7p Cr + t*y for any 7 <= V peT and t £ [0, 1], we deduce 
that 7p er minimizes the following linear functional 

7 G -> ^3 J Tr i|(r) ((ff p ° er )e7 S ) 
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where H® is the mean- field operator defined in (|2.5[) . We subtract the chemical 
potential ep defined above and introduce the functional 

7 € Pper -> F(7) ~ 7A3 Tr ^(r) («* " ^)«7e) 

Notice that since j^s J* r . Tr L 2( F ) (7^)^ = Z for any 7 e "P^ er , then 7p Cr also 
minimizes F on Vp eI - 

For any £ 6 V* , we can find orthonormal functions efc(£, •) € if (r) such that 

(A.2) (H^ = J2^)\e k (t-))(e k ^-)l 

k>l 

each function (£,x) <— > efc(£, a;) being measurable on T* x T. Let us now define 
7° € V pcl by 

( 7 °){(z,y) = £4(0e^;Ml^), 4(0 = { I J[^|j ; ^ 

Saying differently 7 = X(-oo,e F ] (Hpa)- Notice ef was chosen to ensure 7 £ Vp el - 
We now prove that 7 is the unique minimizer of the function F defined above, 
on the set P per without a charge constraint. Since 7 £ Vp el , this will prove that 
7per = 7° an d that 7p er is the unique minimizer of F on "P per . We write 



F( 7 )-F( 7 °) = (2^)- 3 ^Tr i , (r) (( J ff p cr - e ^( 7 -7 ) ? )^ 

= 5>^r 3 / ( A fc(?)-^)((7«e fc (o, efe (e)> c -4(e))^, 



where (•, is the usual inner product of L^(T). Since < 7 < 1 in i 2 (M 3 ), we 
have that < 7^ < 1 on Li(T) and thus (7fefc(£, •), efc(£, •)) G [0,1], for almost 
every £ S T*. Hence, using the definition of 8k(£), 

F( 1 )-F( 1 °)=Y,(^r 3 [ |Afc(0-£F| x |<7eCfc(€,0,efc(C,-)>-«fc(O|dC>O. 
fc>i Jr * 

This shows that 7 minimizes F 1 on Vp el . If now ^(7) = F(7°), then necessarily 
(7£ e fc(£j ")> e fc(£)i ■) = 4(C) I0r almost every £ S T* and any fc > 1, the set {£ £ 
T* I 3k, Xk(£.) — £_f} having a Lebesgue measure equal to zero by [3H Lemma 
2], Using now that the operators 7^ and (1 — 7)5 are nonnegative, we infer that 
7£ e fc(£) — 8k{0 e k(0 for all k > 1 and almost all ^gf. Hence 7 = 7° and 7 is the 
unique minimizer of F. In particular 7p Cr = 7 , i.e. 7p 0r solves the self-consistent 
equation (|2.6|) . 

Consider now another minimizer 7 of the energy £p 0r on Vp CI , we recall that 
p 7 = p 7 o as was shown in [5]. Hence the operators -ffp er and 7 defined above do 
not depend on the chosen minimizer. The above argument applied to 7 shows that 
7 = 7° = 7por> i-e. 7pcr is unique. □ 
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